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Abstract. The quantum dynamical evolution of atomic and molecular aggregates, 
from their compact to their fragmented states, is parametrized by a single collective 
radial parameter. Treating all the remaining particle coordinates in d dimensions 
democratically, as a set of angles orthogonal to this collective radius or by equivalent 
variables, by-passes all independent-particle approximations. The invariance of the 
total kinetic energy under arbitrary d-dimensional transformations which preserve 
the radial parameter gives rise to novel quantum numbers and ladder operators 
interconnecting its eigenstates at each value of the radial parameter. 

We develop the systematics and technology of this approach, introducing the 
relevant mathematics tutorially, by analogy to the familiar theory of angular 
momentum in three dimensions. The angular basis functions so obtained are treated 
in a manifestly coordinate-free manner, thus serving as a flexible generalized basis for 
carrying out detailed studies of wave function evolution in multi-particle systems. 
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1. Introduction 

The challenge for atomic and molecular theory to deal with progressively larger 
aggregates of electrons and nuclei suggests treating them in terms of global parameters, 
in contrast to the usual independent-particle approach Jj], ^, |3|, f|. Globally one may 
represent each configuration of iV particles by a single vector R identified by 3(iV — 1) 
internal coordinates with origin at its centre of mass. -R's modulus R thus represents 
the aggregate's overall size, while its direction R specifies its geometry, i.e., the layout 
of the constituents' relative locations and orientations. Quantum dynamics controls 
then wave functions of R that emphasize -R's evolution from a compact to a fragmented 
structure as R grows. 

In the centre-of-mass frame, the reduced two-body Coulomb problem separates 
in spherical coordinates. Extending the system by adding more particles, the kinetic 
energy retains its three-dimensional spherical symmetry as a subset of its symmetry 
in higher dimensions. This remark suggests concentrating the search for appropriate 
coordinates on the kinetic energy. Hyperspherical coordinates prove suitable for our 
task, by combining the symmetries of particles' kinetic energies into a unified internal 
kinetic energy of the whole aggregate. This procedure ensures an accurate description 
in the compact limit where the aggregate's kinetic energy predominates. 

In a first step, hyperspherical coordinates separate -R's modulus R from its direction 
-R, the latter being represented by parameters analogous to the polar coordinates of 
physical space. The polar-coordinate symmetries, embodied in the familiar formalism 
of angular momentum theory, extend thus automatically to R's (3iV — 3)-dimensional 
treatment. Our proposed perspective of hyperspherical coordinates as extending the 
three-dimensional spherical symmetry to higher dimensions affords advantages far 
beyond the technical aspect of providing a convenient coordinate system. Separating 
the single hyper-radius R from a large number of hyper-angles, and focusing on the 
symmetries under transformations of these hyper-angles, allows one to handle most 
of the many degrees of freedom in a multi-particle system analytically. For example, 
explicit calculation of angular integrals may be avoided in much the same way as 
in ordinary three-dimensional angular momentum theory by replacing such overlaps 
with coupling coefficients derived directly from the symmetry under (ci-dimensional) 
rotations. A second example employs the coordinate-independent representation of ci- 
dimensional transformations || ||, 0, § to construct complete sets of hyperspherical 
harmonics, higher-dimensional analogues of spherical harmonics, without solving partial 
differential equations in any specific coordinate system. This application proceeds 
much as in three dimensions, where all harmonics may be built recursively from a 
single function once the appropriate "laddering" operators are identified. The resulting 
flexibility in choosing centre-of-mass coordinate frames is crucial in the case of a multi- 
particle system with its evolving structure, because no single reference frame proves 
appropriate throughout the entire evolution. 

The hyper-radius R serves as the "evolution parameter" of wave functions 
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^(R; R), whose ^-dependent features evolve with increasing R toward their alternative 
fragmentation terminals. In more detail, the hyper-radius R of an iV-body aggregate of 



masses mi,m 2 , 
given by 



, 771 jV; located at coordinates r±, ?*2, . . . , from the centre of mass is 



R 




N 



M = Y,M l} (1.1) 
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In the special case of charged, point-like constituents (electrons and nuclei) relevant to 
atomic and molecular applications, the Schrodinger equation governing the aggregate's 
evolution takes the hydrogen-like form, in atomic units (a.u.) [|TJ, 



2M 



dR 2 R 2 



z(R) 

R 



R& N -Wy(R; R)=E R^ N - 4 ^(R; R). 



;i-2) 



The factors _R( 3Ar_4 )/ 2 ; proportional to the square root of a hyper-sphere's area with 
radius R, correspond to the familiar factor r of wave functions in physical space. Note 
how (3N — 4)/2 reduces to unity for the H atom, where N = 2. Separation of these 
factors affords non-zero values of \I/ at R = and removes first derivatives from ( |1.2|) . 

The general theory of hyperspherical coordinates (originally conceived for three- 
body scattering problems) dates to the fifties || [10], [11], [12|. Since then, hyperspherical 



methods have been fruitfully applied to a wide variety of many-body phenomena ranging 
from quantum chemistry to particle physics, as illustrated e.g. by several contributions 
in 
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Delves [0, [H| pioneered the method's application to shell-model calculations of 
nuclei. To this day, the hyperspherical approach remains a standard tool in nuclear 
physics, notably in the study of halo nuclei |16|], three-nucleon systems |I7|], as well 
as large nuclei [T^| and even sub- nuclear (quark) structure (For reviews on 



hyperspherical methods in the context of nuclear physics, see e.g. [^D| [H], |22|j.) 

In the seemingly very different context of reactive scattering in quantum chemistry, 
the participating atoms' coordinates have typically been cast in hyperspherical form, 
their motions being governed by effective potentials |23|, |24|, ^5[. Besides molecular 
reactions [|^, [27], [28|, |29|, ^(| [U], ^] the hyperspherical approach applies also to molecular 
and quantum phase effects in chemical reactions 



structure |33J and quantum phase effects in chemical reactions [|34|, p5|| . 

These applications required more general and mathematical investigations on the 



structure and properties of hyperspherical functions |36 



8], |39[ , with a special focus 



on additional quantum-physical features such as anti-symmetrization in fermion systems 



4l| , fl2|1 and internal structures of an iV-particle system in three-dimensional space 



Applications of hyperspherical methods to doubly excited states of two-electron 
systems in atomic physics Jl], resulted in a complete classification scheme for these 
states [0, [15], ^B, K7], based on systematic investigations of the wave functions' structure 
49, 0, 51 1 . Connections with molecular structure [52, |53|] arising from these 



classifications have more recently afforded extending the hyperspherical approach to 
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systems with several heavy particles besides electrons (i.e., to molecules) [Q. In 
three-body Coulomb problems, hyperspherical methods have become computationally 
competitive through the "hyperspherical close-coupling method" ]55 , and more recently 
in the form of "diabatic-by-sector" numerical methods pEfl . Another extension 



employing over-complete basis sets afforded determining accurately resonances of very 
high-lying doubly excited states close to the threshold for double ionization |S7], |55|| . 
Thus, the hyperspherical approach covers essentially the entire energy range from the 
ground state, through the "Wannier-region" around threshold for full disintegration, 



to energies high above this regime (see e.g. [59fl). Besides the two-electron atom and 
generic three-body Coulomb systems, studies have focused on doubly excited many- 
electron atoms pOl Bill, with several forays into atoms with three [B^, |63L |64], BBL B(| 



or even more excited electrons |67|, gg, p9fl . Moreover, the treatment of highly excited 



atoms in external fields |70| , [71|] has extended the list of successful applications of the 
hyperspherical method in yet another direction. Currently, physically adapted Sturmian 
basis sets promise further advances in broader contexts |72) . Reviews on various aspects 



of the hyperspherical approach in atomic physics may be found, e.g., in |73], [74|, [fj|. See 
also JFB) for a discussion of the reliability of the hyperspherical adiabatic method. 

Our equation ( |1.2|) implies going a step further, extending hyperspherical 
coordinates to all constituents, electrons and nuclei, governed by their Coulomb 
interactions. Equation ( |1.2| ) thus represents an exact Schrodinger equation, and all 
calculations proceeding from ( |1.2|) will be completely ab initio [[!], |j, [77], [7^, [79fl . 

We intend in this essay to encompass all types of applications of hyperspherical 
coordinates. Thus the definition fll.lp need only extend over the dynamically relevant 
variables for a given problem, ignoring, for instance, the coordinates of electrons 
belonging to closed shells. In this context note already the interplay between "motions" 
on widely different (time) scales: At each fixed hyper-radius R, a "geometrical" structure 
emerges resulting from the faster motion in the hyper-angular coordinates; the emerging 
structure then evolves on a different scale as R increases. This theme — central to 
the hyperspherical method — serves as a guideline throughout the present paper. As 
a qualitative illustration, consider the following hyperspherical description of a water 
molecule H 2 0, a system consisting of ten electrons, two protons, and an oxygen nucleus. 
In a preliminary step towards ab initio construction of this molecule, a hyperspherical 
procedure would fill first the closed inner shells of the constituent atomic cores, in this 
case only the K shell of 6+ . Recognizing the vast difference between electronic and 
nuclear motion, the electronic motion is "parametrized" by the cores' arrangement. For 
each core arrangement, the electronic motion is analyzed to find the most favourable 
electron distribution. This analysis proceeds in a hyperspherical representation of the 
valence electrons as a single entity consisting at the outset of a group of six electrons 
joined by two single electrons, as suggested by the core charges. Analogous procedures 
for the construction and transformation of such "Jacobi trees" will be outlined in 
section |2|. The particular arrangement of the cores should manifest itself in the hyper- 
radius R e associated with the set of valence electrons alone, once the latter are attached 
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to the cores. The angular coordinates specifying the geometrical arrangement of the 
valence electrons depend parametrically on their own hyper-radius R e , which is in turn 
parametrized by the "size" and "shape" of the nuclei (or cores). The latter is itself 
described by a set of hyperspherical coordinates, {Rn, chjv; @n}, to be introduced in 
section 0, thereby completing the hierarchy of geometrical structures governed by the 
interplay of hyper-radial and angular motions, and of electronic and atomic motions. 
The dynamical evolution of the "shape" characterizing the cores' arrangement should 
minimize the energy of the whole molecule for the correct molecular geometry, while the 
hyper-angles characterizing the distribution of valence electrons should indicate that 
four of them are essentially attached to the 6+ core, with the remaining two pairs 
forming the bonds holding the molecule together. Without attempting to carry out 
the computational details of this description, the present paper introduces techniques 
required for its implementation. Namely, we intend to focus on the universal (i.e. 
coordinate- independent) aspects of hyperspherical coordinates and harmonics, to be 
implemented upon identifying the suitable Jacobi coordinates of a specific system. 

Returning to equation (|1.2| ), we note that features of \& are discrete, owing to the 
finite extent of hyperspherical surfaces, being accordingly represented by appropriate 
quantum numbers and nodal structures rather than by coordinates, much as they are in 
three dimensions; their interpretation will, however, require not only adequate mastery 
of high-dimensional geometry, a subject of this paper, but also of the dynamical interplay 
between light electrons and heavy nuclei. Equation (|1.2|) should serve to calculate 
energy eigenf unctions for any atom, molecule or analogous aggregate. Its solution for 
"collision complexes" formed by colliding molecules will provide the relevant scattering 
matrix directly, as indicated in section ||, by-passing the calculation and study of energy 
surfaces. 



The structure of (|1.2 ) parallels that of the atomic H equation, being actually its 
extension to multi-particle systems. Its first term represents the kinetic energy of the 
hyper-radial motions, its second term that of the hyper-angular motions, and its last 
term the potential energy, which — being the sum of Coulomb interactions among all of 
the aggregate's particle pairs — scales as R~ l . The evolution of ^'s angular part as a 
function of R, foreign to hydrogen, stems from the non-zero value of the commutator 

Equation ( |1.2|) preserves instead hydrogen's invariance under coordinate rotations 
by securing invariance of the kinetic energy operator A R /2Ai under rotations of -R about 
the centre of mass, by mass- weighting the coordinates as described in the following. Sets 
of mutually independent eigenf unctions of A^'s, analogues of spherical harmonics called 
"hyperspherical harmonics", serve to expand ( |1.2| ) into a system of coupled ordinary 
differential equations in the variable R, analogous to those of atomic physics. The 
"effective atomic number" operator Z(R) turns then into a matrix with rows and 
columns labelled by hyperspherical quantum numbers. 

Simple examples of such harmonics in d dimensions have been formulated |5Dl [5T[] , 
the corresponding eigenvalues of A^ having long been known. These examples, to be 
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discussed in section |3|, provide solutions for the second-order eigenvalue problem of A^ 
in specific coordinate systems for d dimensions. Developing an appropriate systematics 
of multi-variable hyperspherical harmonics suited to each system (and thus coordinate- 
independent) constitutes, however, a major objective, to be approached in section [4] on 
the basis of symmetries alone. 

The importance of such a description cannot be overemphasized. The distribution 
of R s multi-dimensional direction Ft should represent specific features of each system 
flexibly, near and far from its centre of mass as well as at intermediate ranges. At 
short ranges, where the system is compact, this distribution should minimize the 
centrifugal effect represented by the eigenvalue of A^ in (|1.2Q . At large ranges, where the 
system fragments, FVs distribution should represent alternative, mutually orthogonal, 
fragmentation channels. Both of these contrasting representations can be achieved in 
terms of flexible "Jacobi (or centre-of-mass) coordinate" sets, interrelated by algebraic 
transformations, the subject of section [2~2"| . 

The present paper develops implications of the quantum discovery that pairs 
of variables conjugate in Hamiltonian dynamics are actually related by Fourier 
transformations, whereby dynamics reduces to kinematics, more generally to geometry 
(whence stems this paper's title). Expanding equation ( |1.2|) in hyperspherical harmonics 
reflects but one aspect of these implications; Jacobi coordinates provide a second aspect. 

The angular Laplacian's eigenvalue in d dimensions reads A(A + d — 2), with 
integer A. Thus, the system of coupled equations resulting from expansion of ( |1.2j ) 
into hyperspherical harmonics parametrized (in part) by the hyper-angular momentum 
A is formally infinite, owing to the infinite range of A, and seems accordingly impractical. 
Key circumstances reduce, however, its size generally to a modest level: 

a) For low- lying channels, the centrifugal term of (|1.2j) , whose numerator rises as A 2 , 
quenches the amplitude of \l/'s components with large values of A to negligible levels 
at small hyper-radius R. 

b) Correspondingly, its Coulomb term (prevailing at large R) has eigenvalues similarly 
spread over many orders of magnitude, its lowest one approximating the lowest 
dissociation threshold for molecules or ionization threshold for single atoms and its 
highest one approaching the threshold for full disintegration of the system. 

c) The range of A values of practical relevance depends thus critically on the energy 
range relevant to each step of evolution. So do accordingly the dimensions of the 
corresponding set of hyperspherical harmonics and of the relevant Z(R) matrix. 



Early calculations [[77], [7^, 79 1 have accordingly shown the range of A values relevant 
at each value of R to be modest, thus affording ready numerical integration of ( |1.2| ). 
Each infinitesimal u dR" step of this integration generates an infinitesimal rotation of the 
^(i?; Ft) wave function by [Z(Ft)/R]dR in the Ft space. The broad range of applications 
envisaged in the present article rests on the power of its underlying recursive procedures 
and on the characteristic aptitude of computer technologies to apply such procedures 
step by step. 
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The vector R, representing the structure of a multi-particle aggregate and 
constructed by recursive procedures, has been resolved above into its magnitude 
R and direction R, corresponding to the aggregate's size and shape, respectively. 
The aggregate's shape, in turn, needs articulating into appropriate parameters— 
multipole moments, for instance — representing structural features of each system, 
a task presenting a challenge to be approached in section |^ in terms of equation 
( |1.2|) 's eigenchannel solutions $ P (^)(-R). An elementary example of such developments 
is afforded by noting that an aggregate's fragmentation elongates its shape, thus 
minimizing its moment of inertia about a symmetry axis. Analogous features should be 
systematically identified and utilized in specific applications. 

The following sections should introduce the reader to analytical tools serving 
to treat atoms and molecules of increasing size: (i) Coordinate systems whose 
dimensionality extends recursively and flexibly, adaptable to particle sets with different 
masses arising in molecular structure and collisions (section 0); (ii) Prototype examples 
of hyperspherical harmonics suitable for multi-electron atomic systems with a single 
heavy centre (section |]); (iii) Systematics of harmonics labelled by eigenvalues of 
commuting operator sets {Hi}, adapted to the evolving structures of atoms and 
molecules at increasing R (section |j). The final section |5] will outline an analytical 
procedure to integrate ( |1.2| ), displaying the evolution of relevant wave functions and 
thereby casting the results of previous algebraic developments |7|, || into a more explicit 
geometrical framework. 

The wide range of symmetry applications, resulting from invariance under 
coordinate rotations and thus relevant for our application to quantum mechanics, is 
currently covered appropriately by |82"| . Terminology drawn from |7|, [8] will appear in 
". . . " , but familiarity with these references is not assumed. 

2. Linear coordinate transformations in higher dimensions 

Physical-space expansions into spherical harmonics hinge on their geometrical and 
kinematical behaviour under coordinate rotations, labelled as the "angular-momentum 
theory". Their extension to larger coordinate sets affords flexibility for displaying 
geometrical and kinematical features of multi-particle systems, at the price of extending 
and elaborating each system's treatment. 

Whereas, in a three-dimensional prototype, rotations about the x-axis, i.e. in the 
2/z-plane, are not independent of rotations about the z-axis because of involving the z- 
axis itself, in a mult i- dimensional setting rotations about different axes are independent 
insofar as they operate in separate planes. Each elementary coordinate rotation is then 
properly identified, in multi-dimensional settings, as occurring within (or parallel to) 
the plane through a given pair of coordinate axes, rather than as preserving a single 
invariant axis. Independent rotations thus occur in non- crossing planes, rather than 
within a single one, e.g., in the xy- and zt-planes in four dimensions. In d dimensions 
the number of independent rotations is readily seen to equal the largest integer not 
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exceeding d/2, i.e., d/2 for even values of d and (d — l)/2 for odd d. This number, called 
the rank of each transformation group 0, HI) is usually indicated as £ = [d/2]. 
Independent Hermitian infinitesimal-rotation operators, corresponding to 



in three dimensions, are indicated generically by Hi, here and in the following. Their 



group's structure and on coordinate choices to be described later. The elementary 
example of the helium atom, consisting of three particles and described in terms of 
six coordinates {xi, y%, z\, x 2 , y 2 , z 2 } with origin at its centre of mass, involves three 
independent rotations represented, e.g., by 



The first two of these expressions are plainly analogues of single-particle's l z \ the 
interpretation of the last one — intermixing two particles' coordinates — remains obscure 
at this point. Note also that finite rotations are familiarly represented by exponential 
functions of infinitesimal operators, as in the example of rotation by an angle <p about 
a z-axis, represented by 



perform in multi-dimensional settings the function performed by l z for three-dimensional 
rotations; their eigenvalues are also integers of either sign, often called "weights". 
After separating the centre-of-mass motion of an iV-particle aggregate, its 3N — 3 
internal coordinates are partitioned typically in three-dimensional sub-sets of single- 
particle coordinates. Sub-sets of independent infinitesimal rotation operators are then 
conjugate to angles ipi = ta,n~ 1 (yi/xi) ranging from to 2tt, and thus effectively 
boundless. These three-dimensional coordinate sub-sets may also include angles &i = 
tan _1 ((a;f + yf) 1 ^ 2 /Zi), 0j's range being restricted by centrifugal potentials near each 
of its poles. Angular coordinates representing tan -1 (ri/rj) of single-particle distances 
from a centre of mass are similarly confined by boundary conditions. For purpose 
of illustration, replace the Cartesian coordinates of helium with six (frequently used) 
hyperspherical coordinates, namely, 




(2.1) 



analytic (differential or algebraic) expression, often analogous to (|2~T|), depends on a 




(2.2) 






1 = 



1,2. 



(2.4) 
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These definitions map the indistinguishability of electrons onto a novel symmetry under 
the reflection a — > 71/ 2 — a. The operators (|2.2| ) take now the form, analogous to the 
last expression in fl2.1| ), 

ctyV d(f 2 ' 9 tan" 1 (cos 6*2 sin a/ cos 6 1 ! cos a) 

Note that sets of commuting operators {-f^i} are subject to coordinate 
transformations among equivalent sets. They are, in fact, suited to represent invariants 
of relevant particle sub-systems. Each of these sets in d dimensions is complemented by 
a much larger set of (generally) non-commuting operators, analogues of l x ± il y in three 
dimensions that raise or lower the m quantum number of spherical harmonics Yi m (8, ip), 
respectively. Thereby one reaches the total of d(d — l)/2 (increasing quadratically 
with d) unrestricted linear infinitesimal transformations in d dimensions. Most of these 
operators characteristically involve coordinates of different particles, thus correlating 
their motions. E.g., combinations of equation ( |2.1| )'s analogues 

would raise and/or lower eigenvalues of operator pairs {Hi,Hj}. [The upper indices 
of the J symbol denote the physical-space components, while its lower labels signify 
particles' indices, together specifying the relevant variable pair in the <i-dimensional 
space.] 

This and analogous considerations enlarge the scope of our study considerably, 
yielding a total of d(d — 1)/2 = (3iV — 3)(3iV— 4)/2 infinitesimal operators for N particles 
in d = 3(N — 1) dimensions. The resulting wealth of operators will be introduced here 
and developed later. As each Hi involves two coordinates, odd values of d imply that 
one coordinate, often labelled by 0, fails to be included in any of the I H^s, even though 
contributing to the set (pT6|) . 

To establish contact with the relevant mathematical theory of Lie groups and Lie 
algebras 0, [8], 82], we observe that the set of linear infinitesimal rotation operators 
in d dimensions considered in this section forms the "special orthogonal group" of 
transformations, designated as SO(d); the term "special" referring to exclusion of 
dilations. Quantum mechanics extends this group to include the analogous complex 
transformations forming the unitary group SU(d) of transformations that preserve the 
complex Hermitian form XkX* k instead of its real quadratic analogue. The full 
unitary group U(d) results from adding multiplication with a complex phase exp(ia) 
to its "special" version SU(d). The real part of the Hermitian form is preserved by 
orthogonal transformations, its imaginary part by "symplectic" transformations that 
preserve bilinear forms odd under permutations of their elements, such as the spin 
invariant ^i^Ix^ — u l/2 u -i/2- 
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2.1. Symmetry under rotation reversal; ladder operators 

Reversal of a rotation's direction is represented, e.g. in fl2.ip, by switching the sign of 
the coordinate y or of the imaginary unit, thereby reversing the sign of the operators l z 
and of their eigenvalues m. Symmetry under this reversal has been high-lighted in |7], 
by replacing pairs of coordinate labels (xi,yi) with pairs (xi,X-i), thus replacing yi in 
( p.6|) . The index i runs thus over (1, 2 . . . ,£), being complemented by an x$ for spaces 
of odd dimension d. Spherical coordinates include then £ azimuths ipi, with odd parity 
under rotation reversal, and i polar angles 8i, changed by this reversal into tt — 

Here we preserve the familiar notation with coordinate pairs pointing 
out that the rotation reversal is often complemented with the reflection through 
the coordinate plane xz, which automatically reverses the sign of The resulting 
combination reverses the handedness ( "chirality" ) of each particle's space coordinates 
(x,y,z). 

Note first that, whereas the physical space operators (l x , l y ) change by unity the 
eigenvalues m of l z , each of the iV-particle operators ( |2.6| ) shifts the eigenvalues {mi, mj} 
of an operator pair {Hi,Hj}. Whereas combinations (l x ±il y ) act as "ladder" operators 
raising or lowering the eigenvalues m of l z by unity, combinations of four among these 
iV-particle analogues ( |2.(j| ) raise or lower eigenvalue pairs simultaneously. The resulting 
rather elaborate classification of operators became a central feature of the algebraic 
treatment [[I], [J; we shall follow a more direct approach. 

Recall how the ladder operators of physical space, l x ±il y , emerge as non-Hermitian 
combinations of the Hermitian pair (l x ,l y ). Similarly, non-Hermitian "raising" and 
"lowering" conjugate operators, designated here generically as (aJ,a), respectively, 
may be viewed as combinations of two pairs of Hermitian operators, a' + a and 
i(a^ — a), symmetric and antisymmetric, respectively. Recall also that the physical- 
space Hermitian operators (l x , l y ) are antisymmetric and symmetric, respectively, under 
reflection through the xz-plane, according to standard "Condon-Shortley" conventions. 
A heuristic approach to constructing ladder operators might thus start by identifying 
combinations of operators ( |2.6|) that are symmetric and antisymmetric, without 
resorting to Lie-algebra procedures. 

Consider now the example of infinitesimal operators pertaining to the four- 
coordinate set {xi,yi,Xj,iij}. Besides the commuting pair {Hi,Hj}, denoted here as 
{(xi, yi), (xj, yj)}, the set of operators includes the ones from (|2.6| ), labelled (xi,yj) 
and (xi,Xj), respectively, as well as the additional two (xj,yi), (yi,yj), totaling six 
operators, four of which intermix i and j coordinates. This last sub-set gives rise to 
two pairs of symmetric and antisymmetric Hermitian operators, and thence to two 
pairs of desired non-Hermitian operators, one of them "raising-raising", designated as 

(++), together with its conjugate ( ), and "raising- lowering" operators (H — ) and 

( — h). Commutators among these operators generate other operators of the set, much 



as commutators among the I components do, as detailed in [Appendix A| . This particular 



example, with ladder operators appearing simply as combinations of two components 
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acting on and rnj in the way familiar from SO (3), is not generic for SO(d). Instead 
its structure is due to the well-known feature of orthogonal groups that 50(4) factors 
into a pair of SO (3) sub-groups. A related example familiar in physics is afforded by the 
factoring of the (proper) Lorentz group S0(3, 1) into 5X7(2) x SU(2), with "spinors" of 
opposite chirality [^, Chapter 11]. Reference J7| derives these results algebraically. 

Returning once again to the physical-space operators l x ± il y , recall how they 
transform elements of a spherical harmonics set {li m } into one another. They serve 
further to identify a set's range, e.g., by causing its "highest weight" element Yu(9, ip) 
to vanish when acted upon by l x + il y . The ladder operators outlined above will 
perform an analogous role for the much larger and multi-faceted sets of hyperspherical 
harmonics. To this end one may combine the operators ( p.6|) first into pairs symmetric 
and antisymmetric in their indices, 

= ^ M? ± ^0 = ^ W T J??) , (2.7a) 

since Jjf = —Jtf- These Hermitian operators are then paired into non-Hermitian 
operators, 

The heuristically introduced operators of the type ( |2.76|) fail to act correctly as ladder 
operators. Nevertheless, examples of ladder operators with a similar structure arising 
as superpositions of four operators (|2.6| ) will appear in ( 4.11o| -c), playing a key role in 
section [|. [Combinations (|2.6|) that include instead the unpaired variable x of odd- 
dimensional systems do not lend themselves to the symmetrization ( |2.7q| ), being thus 
more nearly analogous to the three-dimensional (l x ,l y ).] 

The numbers of infinitesimal operators ( |2.6| ) and of the resulting ladder operators 
increase quadratically with the numbers of particles and of the corresponding 
coordinates. It will turn out in following sections, however, that a number I of linearly 
independent "raising-lowering" operator pairs, equal to the number of commuting 
operators ifj, suffices to generate complete orthogonal sets of hyperspherical harmonics. 
Each of those sets corresponds to a choice of relevant coordinates and of the {Hi} set. 

The coordinates and their infinitesimal rotations, developed thus far in this section, 
would intermix in their dynamical applications with the inertial effects of the mass- 
differences among various particles. These complications can, however, be removed 
by appropriate mass-weighting of the coordinates as anticipated in section [l] and 
implemented next. 



2.2. Jacobi coordinates 

Mass-weighting of coordinates has served in section [l] to define the hyper-radius R 
conveniently, contrasting it with the angular coordinates represented by -R. Analogous 
devices serve to weight appropriately the components of R pertaining to particles with 
different masses Mi, making them homogeneous, and thus removing the Laplacian A^'s 
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explicit dependence on single-particle masses. To this end, generic sets of "Jacobi 
coordinates" have been introduced long ago flT2"|| , replacing the single-particle r, by 
vectors £ 4 with the mass-weighted dimension mass 1 / 2 -length. 

Alternative sets of Jacobi coordinates occur, reflecting alternative groupings of 
particles, properly weighted by the mass of each group through linear transformations 
with dimensionless coefficients. These sets are in turn interconnected by dimensionless 
linear transformations, each of whose steps amounts to rotation in one plane. Handling 
of Jacobi coordinates becomes thus laborious even though each step be elementary. 

Jacobi coordinates prove essential by identifying each fragmentation channel 
through a particular Jacobi tree separating at its base into two branches corresponding 
to the relevant fragments. One thus displays the evolution of each particle-aggregate 
toward a specific fragmentation channel by the structure of relevant tree-shaped Jacobi 
coordinates. The label "Jacobi tree" reflects the evolution of a multi-particle system 
fragmenting (i.e., "branching out") from a trunk into separate systems. [Developing a 
multi-particle wave function toward one among its alternative fragmentations presents 
instead an upside-down view of that tree.] 

We show below a few simple prototype Jacobi trees, whose upper endings 
correspond to single-particle labels. Permutation of two particles is represented by 
rotating (a) about its trunk. Changing of (b) into (c), often referred to as transplanting 
of branch 2, corresponds to the prototype transformation of hyperspherical harmonics. 
Generic transformations resolve into sequences of permutations and transplantations. 
Tree (e) is often called "canonical". Multi-particle systems are represented by 
correspondingly articulated trees. 

1 2 3 4 5 

(e) 

The Jacobi trees just introduced serve to characterize multi-particle systems by 
their hierarchy of composition, i.e., by indicating the order in which particles are joined 
to form sub-complexes of the entire aggregate. The labels thus refer to particle indices. 
A similar concept re- appears in a further context: Analogous trees illustrate coordinate 
systems and their transformations, their branches labelling appropriate angles in FVs 
decomposition |2D|. Rotation of tree (a) by an angle < tp < 2tt represents then 

a simple rotation about an axis. Variation of a coordinate < 9 < it maps onto the 
angle between two branches. 

Alternative sets of Jacobi coordinates R correspond thus to different tree structures. 
This circumstance adds further elaboration to our procedure, yet serving to display 
whole aggregates' evolutions. These aspects have not been apparent in the initial 
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applications of the present approach, dealing with very few particles, even though 
necessarily underlying the treatment of any multi-particle system. Each transformation 
of hyperspherical harmonics resolves accordingly into a transformation from one to 
another set of Jacobi coordinates and a transformation of the corresponding harmonics. 

Note the "hierarchical" aspect of Jacobi-tree construction, which adds particles 
sequentially, contrasting with the "democratic" view of the multi-particle coordinates 
leading to the quadratic increase of the number of fl2.6|) operators as a function of the 
particle number N. 

Once again: Each step dealing with Jacobi coordinates is elementary, but the 
number and combinations of different steps are large, a characteristic generally 
encountered in computer operations requiring adequate strategy and planning. Such 
operations may properly articulate into successive phases. We anticipate, for example 
in dealing with molecules, to build first each atom's inner shells independently, by 



Cavagnero's procedure |67], |69|, combining later the resulting atomic ions with the 
residual atomic electrons. 

Similarly, the following sections present first alternative combinations of three 
particles, lying in a plane with their centre of mass, into alternative pairs of mutually 
independent collective coordinates followed by transformations among these pairs. 
Combinations of larger particle sets into further £'s will be dealt with next, utilizing 
analogous procedures recursively. Even more extensive procedures will hinge on 
experience in treating large multi-particle aggregates. 

2.2.1. A three-particle prototype The positions r*, i = 1,2,3, of a three-particle set 
identify a plane where their centre of mass also lies. These positions are represented 
by these co-planar vectors, but their internal kinematics involves only two independent 
Jacobi coordinate vectors, £, with 3(iV — 1) = 6 degrees of freedom. 

The standard procedure for constructing Jacobi coordinate vectors considers first 
the positions of two among these particles, 1 and 2 here, with masses M\ and M 2 and 
positions 7*1 and r 2 . These input data combine into a first Jacobi vector weighted by 
the square root of the pair's "reduced mass" Myi = M\Mij{M\ + M 2 ), namely, 



£12 = y/M^in - r 2 ), (2.8) 

whose centre of mass lies in the plane at 
Mxn + M 2 r 2 

ri2 = M 1 + M 2 • (2 ' 9) 
The next step combines the first pair, with mass Mi + M 2 and centre of mass 
position 7*12, with the third particle lying at r%. This step is performed in accordance 
with (|2.8|) yielding the second Jacobi vector, 

^i2,3 = V /M ^(^i2-r 3 ), (2.10) 

with the reduced mass 

(M 1 + M 2 )M 3 

' d M 1 + M 2 + M 3 K J 
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[The comma-separated sub-scripts indicate the particle sub-complexes to be joined.] 

An additional feature relates the Jacobi vector ( |2.10| ) to its alternatives 
corresponding, e.g., to the permutation of indices 1 and 3, i.e., (1,2,3) — > (3,2,1), 
yielding the two-dimensional vector rotation 

{£s2> £32,1} = {£i2 cos /3-£i2,3 sin A £12 sin/ 3 + £12,3 cos/?}, (2.12) 

by the angle 



Analogous kinematic rotations correspond to cyclic permutations of indices. 



2.2.2. Extension to multi-particle aggregates The formulation of equation Q2.10| ), with 
elements from ( |2.8|) , has clearly recursive character. It implies that any pair of Jacobi 
vectors, {£ p , £ q }, representing two sub-aggregates of particles centred at r p and r q with 
masses M p and M q , respectively, combines effectively into a single vector 



*pq 



(r v -r a ). 2.14 

M p + M q p q! y ! 

Similarly, restructuring of any "Jacobi tree" diagram, which represents a specific 
sequence of particle combinations forming an aggregate, resolves into sequences of 
vector-pair {$, p ,£ q } rotations within a plane, analogous to that represented by ( |2.12| ). 
Such restructurings have been discussed amply in []20| under the name of "timber 
transformations" , the word "timber" being suggested by association with "Jacobi tree" . 
The simple underlying principle, stated in that reference, lies in the feasibility to resolve 
any rotation in mult i- dimensional spaces into a sequence of plane rotations, a feature 



familiar for the three-dimensional rotations of physical space. Appendix B| exemplifies 
this procedure. 

Transformations between different Jacobi trees prove highly relevant to our subject 
of atomic and molecular few-particle systems for the following reason. The Coulomb 
coefficient Z(R) results familiarly from contributions proportional to the reciprocal 
distances between the N(N — l)/2 particle pairs. Flexible sets of Jacobi coordinates 
afford treating each of these distances as a single coordinate, to be combined with 
others, thus avoiding the familiar need to expand each term of Z(R) into a multipole 
series. Thereby determining Coulomb interaction matrix elements reduces to calculating 
integrals over 1/r rather than over 1/ \ri — r 



3. Sample hyperspherical harmonics 

The familiar spherical harmonics, Yi m (8,(p), serve as tensorial base sets for (21 + 1)- 
dimensional transformations induced by rotations of the physical-space coordinates. 
In multi-dimensional contexts analogous base sets of hyperspherical harmonics serve 
the same purpose. The name "harmonics" identifies them as eigenfunctions of the 



Geometry and symmetries of multi-particle systems 



15 



angular Laplacian operator A^ in the (d — l)-dimensional space of R with eigenvalues 
-A(A + d-2). 

To understand the term d — 2 in this eigenvalue formula, note first that it reduces 
to unity for d = 3 yielding the familiar eigenvalue 1(1 + 1) of the squared orbital angular 
momentum. The unit in this expression corresponds to the single angular coordinate 9 
that accompanies the angle <p in polar coordinates. Its contribution to the eigenvalue, 
namely I, corresponds to the "zero-point energy" , u hF , of a unit-mass particle oscillating 
along the 9 coordinate in the centrifugal field generated by its rotation along ip with 
quantum number /. The occurrence of d — 1 dimensions for the vector R raises the 
number of its coordinates, besides (p, from unity to d — 2, thus accounting for the 
eigenvalue term d — 2. 

As the pair of angles (9, p) identifies a direction of physical space, an equal number 
of indices (l,m) identifies a harmonic belonging to a (21 + l)-dimensional set, with the 
magnetic quantum number m labelled as a "weight" and I in the role of "highest weight" . 
Extending this parametrization to rf-dimensional spaces requires us to describe sets of 
hyperspherical harmonics Yx^R), where A replaces the "highest weight" I, the vector 
H represents a set of d — 2 complementary labels, and R a corresponding set of d — 1 
angles. 

This extension provides a main tool for the quantum mechanics of multi-particle 
systems, as indicated in section [1] anticipating the relevance of the hyperspherical 
harmonics and of their treatment in |77], [80]]. Recall that the position vector R of 
an iV-particle set (in its centre of mass frame) has d = 3(N — 1) dimensions. The wave 
functions ^(R) of such a system, envisaged in section [TJ, are conveniently expanded in 
hyperspherical harmonics in analogy to expansions in spherical harmonics ||77|| . Their 
Schrodinger equation reduces similarly to a system of coupled ordinary differential 
equations in the "hyper-radius" R. 

The symbol for hyperspherical harmonics, Yxn(R), replaces the index / of spherical 
harmonics by the index A corresponding to the eigenvalue — A(A + d — 2) of the (d — 1)- 
dimensional angular Laplacian A^. This index also stands for the degree of the 
homogeneous "harmonic" polynomial products R^Yx^R). The second index fj, replaces 
the index m of spherical harmonics with a corresponding set of d — 2 parameters (the 
dimension of R less 1) that identify a specific harmonic of degree A. 

Expanding wave functions of a multi-dimensional R utilizes "complete orthogonal 
sets" of hyperspherical harmonics. Completeness is achieved by extending the range 
of A adequately. Only a finite set of harmonics proves, however, relevant at any finite 
value of R, higher values being effectively excluded by a generalized centrifugal potential 
at small R, as noted in section [D. [This potential, involving the A parameter, includes 
contributions from derivatives of variables corresponding to the a coordinate in ( |2.4| ) and 
representing the quantum mechanical resistance of particles to compression by boundary 
conditions.] As R increases more and more hyperspherical harmonics start contributing 
to the relevant wave functions, requiring adequate frame transformations to reflect the 
appropriate fragmentation channels, as outlined in section |5|. Orthogonality requires 
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identifying, for each A value, an adequate set of vectors fj, corresponding to the simple 
set of integer values |m| < I of the spherical harmonics and to the relevant set of R 
components. Higher dimensionality implies here more elaborate sets of \x vectors. 

Generating these sets in d dimensions, labelled by (d — 2)-dimensional vectors 
H, relies in essence on separating the Laplacian's variables. In the three-dimensional 
prototype Yi m (9,ip), the m label arises as an eigenvalue of the l z operator representing 
the number of nodes of the corresponding eigenfunction sin' m ' 9e imip , whereas the 
remaining I — \m\ nodes pertain to the 9 variable. For hyperspherical harmonics, sets of 
section £ commuting operators if, may provide corresponding quantum numbers 
and eigenfunctions. The residual A — J2i l m «l nodes would then pertain to the d — £ — 1 
residual variables, analogues of 9. 

Whereas alternative orthogonal sets of spherical harmonics pertain to alternative 
orientations of the z coordinate axis in three-dimensional space, analogous sets of 
hyperspherical harmonics pertain to alternative selections of £ commuting operators 
{Hi}, not necessarily based on a coordinate set as they were in section |[ Whereas 
spherical harmonics Yi m (9,(p) depend on the "longitude" <p with \m\ "meridian" nodes 
and on the "co-latitude" 9 with I — \m\ "parallel" nodes, sub-dividing their plots 
into separate "lobes" , hyperspherical harmonics depend — for particular coordinates — on 
"longitudes" tpi, each with \rrii\ nodes, and on the nodal distribution in the remaining 
coordinates. Additional components of /x pertain to alternative partitions of A that 
delimit the range of the parameters \rrii\ as well as of additional coordinates. Note 
how the number of \x components increases linearly with the number N of particles, 
contrasting again with the operators Q2.6| ), whose number increases quadratically with 
N. 

Consider now how the features of hyperspherical harmonics bear on equation ( |1.2| )'s 
expansion: The value of in its centrifugal term depends only on the parameter A 
of each harmonic, whereas the coefficient Z(R) of its Coulomb term resolves for an N- 
particle set into N(N — 1)/2 terms ~ R/\ri— Vj\, cast as matrices in the {A^t} basis. The 
flexibility afforded by selecting the £ operators Hi and the remaining d — l — 1 coordinates 
should serve here to avoid, or at least minimize, resorting to multipole expansion of each 
term. 

We anticipate that coordinate rotations in d dimensions transform generally any 
harmonic Y\^(R) into a superposition of the harmonics of its whole orthogonal set 
with coefficients D^, analogous to those that serve to transform spherical harmonics 
and form a "representation" of the rotation group. This feature rests on the multi- 
dimensional rotations' aptitude to resolve into sequences of two-dimensional rotations, 
thereby affording us to express any D^, coefficient in terms of more familiar three- 
dimensional Wigner's Euler-angle functions d®, m . The occurrence of alternative sets 
{Xx/x} with the same A value reflects the multiplicity of the infinitesimal operator basis. 
Systematic classifications of alternative sets {Yx^}, appropriate to the number and 
structure of their coordinates, shall rest on the sub-group chains of the relevant group, 
as outlined below. With this background we describe now two sample harmonics' sets. 
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3.1. A generic structure 

Laplacian equations for hyperspherical harmonics in d dimensions, 

[A A + A(A + d— 2)]Y Xll (R) = 0, (3.1) 

lend themselves to solution by separation of variables, since A^ amounts to a sum 
of terms f{R)(d/dxi)gi{R){d/dxi) with metric coefficients f(R) and g%(R). Their 
prototype example is afforded by the Schrodinger equation for the He atom, with N = 3 
(1 nucleus and 2 electrons), d = 3(N — 1) = 6, and with three polar coordinates 
yielding 



1 9.o 2 d 
An = -^ — —sin a cos a— — h 



sin a cos" 1 a act act cos^ ct 



9 . „ 9 Id 2 
sm ■— — + 



sin 0i 90i d6 1 sin 2 0i 9^ 



sin 2 ct 



19 3 1 9 2 

sm 2 777 — h 



sin 2 90 2 90 2 sin 2 9y? 



I / 9 2 Z 2 Z 2 

' +4 1 r^-)sin2ct, (3.2) 



sin2ct \da 2 cos 2 ct sin 2 ct / 

whose form on the last line arises from renormalizing the harmonics by the volume 
element sin2ct. 

A familiar approach Jl|, ||, [3|, [|] to solving ( j3+l) with the angular Laplacian (|3.2| ) 
assumes first F A/X to depend on (ipi, <p 2 ) through a factor exp(imi<£>i + im 2 {p 2 ), whereby 
each (d/dfi) 2 element of ( |3.2| ) amounts to — m\. Thereafter each of the square brackets 
in ( |3.2|) reduces to —li(k + 1), provided depends on each 0j through the associate 
Legendre function sin m * 0; P z mi (cos0j), with the relevant "highest weight" ij thus limiting 
the range of \mi\. The residual operator on the right of ( |3.2j ) has then eigenvalues 
— A(A+4), with A partitioned as Z1+Z2+271, and with the "Jacobi polynomial" eigenvector 
pfa+ 1 / 2 > ll + 1 / 2 )[ cos 2a), according to (22.6.4) of | 83]| . [The factor 2 multiplying n and ct 



in these expressions stems from the exponents and coefficients in the functions of ct in 
Q3.2D 1. The eigenvector of (|3.1j ) reads thus 



Y XfM (R) = cos' 1 a sin /2 ctPi' 2+i ' 1+ ^ ) (cos2ct)r iimi (0 1 , <px)Y hma (9 2 , <p 2 ), (3.3) 

with the expected five-component n = {n, / 1; m 1; l 2 , m 2 }, and with the spherical 
harmonic factors Yi im .(9i,ipi). 

The set of harmonics Q3.3| ), with a given value of A, consists of a number of elements 
w, the dimensionality of the relevant space. This number depends on the number of 
alternative partitions of A into (h,l 2 ,n) consistent with the relation A = li + l 2 + 2n, 
and on the alternative 21 i + 1 values of each m^. The prototype example of A = 2 leads 
to the partitions: (2, 0, 0), (0, 2, 0), (0, 0, 1), (1, 1, 0) and, in turn, to 5 + 5 + 1 + 9 = 20 
harmonics. Different partitions of A yield hyperspherical harmonics with alternative 
nodal patterns, reflecting alternative sharing of rotational kinetic energy in different 
modes. The algebraic determination of w is discussed in |7|, || [82| and outlined in 



section |[ The general expression || for w reduces in the present case (d = 6) to 
(A + 3)(A + 2) 2 (A+ 1)/12. 
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Notice how the eigenvalue parameters (h, l 2 , 2n) play the role of "weights" in ( |3.3| ), 
each of them amounting to the "highest weight" for a sub-group of the rotation group 
SO(d = 6), thus contributing to the "highest weight" A = h + I2 + 2n of the 5*0(6) 
harmonics ( |3.3| ). The symmetry under the sign reflection (m — > — m) (cf. section |2.1| ) 
manifests itself in these harmonics not only through the symmetry of each factor Yi imi , 
but also through the parity of the polynomial P n , which is even or odd for even or odd 
values of the "pseudo weight" n, combined with interchange of its upper indices. Note 
also how n substitutes for the weight m 3 , an eigenvalue of the third operator Q2.5Q , 
which fails to commute with the operators If even though commuting with (d/dtpi) and 
(d/d V2 ). 

The harmonics (|3.3| ) provide, of course, a basis for expanding correlated electron 
wave functions in helium, in the form ty(R;R) of section [l], with the direction 
R = {a, #1, ipi, 9 2 , f 2 } |77|, IZ11- The expansion coefficients ^ n,h,mi,h,ma{R) represent 



then desired features of the relevant eigenfunction ^{R; R). This representation does, 
however, emphasize single-electron aspects of the state through its parameters 
rather than through the global features anticipated in section [1]. 

Note, on the other hand, that the procedure presented above to construct the 
harmonics ( |3.3|) applies recursively to (N > 3)-particle systems, as should the previous 
comments, with the following qualifications: Extending ( |3.3| ) to iV > 3 particles involves 
N — 1 spherical harmonics Yi imi (9i, <pi), initially independent of one another, and N — 2 
Jacobi polynomials in cos 2a j with the full set of 3N — 4 variables. The index pairs of 
these polynomials, replacing the single-particle (U + 1/2) of (|3.3j), reflect however the 
total angular momenta, L or J, of paired particle sub-sets as determined by relevant 
hierarchical additions of single-particle momenta . Note also how the structure of the 
harmonics p.3j ), and of their higher-dimensional analogues, exploits the commutability 
of the base operators {Hi} and of their corresponding sub-group structure, with a 
notable exception: The angle a as defined in fl2.4Q does not coincide with the arctan of 
the third operator Hi of (|2.5| ). 



3.2. An alternative structure 

The particular sub-group structure of ( |3.3| ) stresses the single-particle features of each 
system, contrary to section emphasis. We turn now accordingly to [8IJ 's alternative 
construction of harmonics that emphasizes different elements, dealing with d = 3N — 6 
unspecified coordinates besides a spherical harmonic Yi m (8,<p), without reference to 
single-particle positions. This approach thus utilizes a single eigenvector e linv of a single 
operator Hi, contrasting with the full set occurring in (|3.3|) and its extensions to iV > 3. 

As a preliminary to complementing a single spherical harmonic Y/ m with additional 
variables, reference |8(| views the Legendre polynomial P/(cos#), invariant under 



rotations about the the particular Gegenbauer polynomial C"(cos#) with 

a = 1/2 (as defined in ||83|| , Table 22.6) extending it from 3 to d = 3(N — 1) dimensions 
by setting its parameter a at d/2 — 1, and replacing I with A (d = 6 in our example). It 
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also replaces cos 9 with the scalar product of unit vectors R- R', where R' corresponds 
to the reference axis z. The hyperspherical harmonic 

C d x /2 ~\R- Ft!) (3.4) 

is thus invariant under d-dimensional rotations of R about a fixed axis R' (i.e., rotations 
labelled by a single parameter if), as well as under rigid rotations of the pair R ■ R! . 
Note that Gegenbauer (as Legendre) polynomials consist of only even or of only odd 
powers for even or odd values of A or I, respectively. 

As the spherical harmonics Yi m are generated from Pi(cos8) by operators l x ± il y , 
harmonics Yx^R) are generated from (|3.4j ) by infinitesimal operator analogues of l x ±il y , 
each of them combining a derivative, which lowers by unity the degree of their operand, 
with a compensating factor linear in a coordinate, thus replacing a nodal line of the 
operand with a different coordinate's nodal line. [Typically, the operator l x + il y , as 
applied to Y\ m with non-negative m, is equivalent to e lv sin 0(5/5 cos 0): Its partial 
derivative reduces V5 m 's polynomial dependence on cos 9 by one degree, thus removing 
one parallel-line node, its first factor adds instead a meridian-line node implied by the 
vanishing of its real or imaginary parts. To trace out this effect more explicitly, recast 
the operator in the form 

, ., .( d 5 \ / 5 5\ . . , 5 f d .5 

l x + il y = -l y— - z — + z- x— = -{x + iy) — + z — + i— 

\ oz ay J \ ox oz J az \ox ay 

The first term on the right of this expression, when applied to the function Yi m (9, if) 
with a non-negative m value, raises its m index by one unit through the combined action 
of its two factors: The first factor, x + iy = r sin 9e lip combines with the same factor 
within Yi m , thus raising by one unit the number of its meridian-line nodes as well as the 
exponent of its sin m 9 factor. Its second factor (d/dz), equivalent here to r~ l (d/dcos9), 
reduces by one unit the degree of Y/ m 's polynomial dependence on cos 9, thus suppressing 
one of its parallel-line nodes. The second term on the right of the operator's expression 
cancels instead the function Y\ m {9, if) by setting it to zero.] 

A direct analogue of the successive action of l x + il y operators on P^cos^) would 
apply to the Gegenbauer polynomial (|3.4|) a sequence of corresponding operators acting 
on components of R. A novel aspect emerges, however, at this point, since R has 
several components orthogonal to R! acted upon by alternative operators in alternative 
sequences. Such alternative sequences will occur in section ^, whereas reference [B(J 



avoids them by treating all R components uniformly, besides their physical-space sub- 
set R = {sin 9 cos if, sin 9 sin if, cos 9}. 

To this end, reference [|80]| follows the frequent practice of complementing an initial 
two-dimensional component (cos ip or sin if) by multiplying it with the factor sin 9 and 
adding a further R component cos 9, whereby R retains its unit magnitude. Iterating 
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this extension d — 3 times yields the canonical set of R components 

COS#i 



sin Q\ cos 6- 



2 



: (3.5) 
sin 6*i sin 6 2 ■ ■ ■ cos Qd-2 

sin 9i sin 9 2 ■ ■ ■ sin #^-2 sin (p 

sin #i sin 6> 2 ■ ■ ■ sin #d_2 cos yj. 

For purposes of orientation consider that, if all R components had comparable 
magnitudes, each value of cos^j would be of order 1/d, and hence each 9i close to 
7r/2. After elimination of RcosOi, the factor Rsm8\ of all successive components would 
represent their total magnitude. The successive factors sin 9i contribute then to reduce 
the effective residual components of R progressively. 

The set of integer components of the vector label /x, non-negative and of decreasing 
magnitude, is similarly indicated by 

/Li = {//i,/U 2 , -Hd-i}, (3-6) 

with fid-i = \ m \-, m being the multiplier in the harmonic's phase imp. The difference 
between successive components (|3.6|), fij — > 0, represents the number of nodes 
in the j-th harmonic's dependence on cosOj, including the nodes implied by the 

phase factor e imip . The set of harmonics compatible with these specifications, for the 
same prototype values A = 2 and d = 6 as for the harmonics (|3.3j ), has likewise 20 
elements. 



The resulting hyperspherical harmonic, equation (3.69) of |p0 |, consists thus mainly 
of products of d — 2 Gegenbauer polynomials, 

d-2 

K AM = N Xtl H(sin Oj^CfitjEl (cos e 3 )e^. (3.7) 
i=i 

The combined degree of this harmonic, equalling its total number of nodes, amounts 
to the X^=i(/% — A*i+i) = \ including the iid-\ = l m l nodes attributable to the e 1 " 1 ^ 
factor. The a index of each Gegenbauer polynomial consists of two terms, ctj and 
fJij+i, the first of which, aj = (d — j — l)/2, corresponds to the actual dimensionality 



of all the ( |3.7|) factors with indices larger than j. The second term, fij+i, represents 



the additional effective dimensionality attributable to the factors (sin#j) Mi+1 [whose 
combination with C^ 3 amounts to an analogue of the associate Legendre function 
P/"(cos0)]. The remaining factor of ( p.7| ), Nx^, represents the normalization factor 
contributing to orthonormalize each set of harmonics ( |3.7[ ) with equal A and alternative 
fi indices. 

Each factor (sin0 3 -) w+1 , which decays rapidly as dj approaches a pole (at or it) 
corresponds to the factor (sin^)' m ' of associate Legendre functions, which reflects the 
centrifugal potential generated by its factor e imip . Indeed the exponent /ij+i equals 
the total number of nodes in ( |3.7|) factors with indices larger than j, which contribute 
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a centrifugal potential to the equation governing C^'.-fi+l', each difference fij — 
corresponds to a separation parameter 1^(1^ + 1) in the construction of the harmonics 

In conclusion, the present construction of hyperspherical harmonics has followed 
two approaches: (i) Solving the (angular) Laplace equation by separation of variables, 
leading to the harmonics ( |3.3| ) and extensible to higher dimensions; (ii) Constructing 
the harmonics (|3.7f ) by a sequence of Gegenbauer polynomial factors (also separating 
variables), complemented by factors (sin analogues of the (sin#) m of spherical 
harmonics and similarly generated by applying infinitesimal operators sin 9(0/ d cos 9) 
to the invariant Gegenbauer harmonic fl3.4|) . 

These approaches differ in their coordinates as well as in their dynamical 
implications: No physical specification of the R components 9j has occurred in the 
approach (ii). Approach (i) has relied on the single-particle coordinates and on their 
interrelations introduced in ( |3.2|) , namely {0 < a < tt/2, < 9i < tt, < (pi < 2ir}. 
[The limited range of a coordinates reflects their definition through ratios of non- 
negative variables.] Approach (ii) has utilized, in (|3.5| ), ratios of R components restricted 
to real values of either sign, in addition to a single complex phase, without explicit 



reference to single-particle coordinates (considered elsewhere in [^0|). The occurrence 
of two (or more) complex-coordinate phases in the approach (i) utilizes the commuting 
operator set {Hi}, a dynamical element foreign to the approach (ii). Intermediate 
approaches, utilizing that set partially, appear readily accessible. 



4. Classification and construction of hyperspherical harmonics 

Section ^| introduced hyperspherical harmonics for two different sets of coordinates, 
relying on single-particle features to a different extent. Both sets of hyperspherical 
harmonics are characterized as "harmonic polynomials", i.e., as eigenfunctions of the 
angular Laplacian, obtained directly by separation of variables. Their construction thus 
ties each of these harmonics to a particular choice of coordinates. 

By contrast, sections [3] and ^| repeatedly stressed the need for flexibility in 
the choice of coordinates to describe the evolution of an atomic or molecular 
complex from its compact to its fragmented states. Having familiarized the reader 
with properties of hyperspherical harmonics in the preceding section, we now 
introduce harmonics essentially frame-independent, thus by-passing extensive frame 
transformations necessary for harmonics tied to a particular set of coordinates. We 
describe here a procedure to generate complete sets of harmonics independent of the 
choice of coordinates. Consequently, the same set of harmonics serves throughout the 
entire evolution process, frame transformations reducing to the task of expressing the 
harmonics in whichever coordinate system appears appropriate at any given stage of the 
evolution, without change to the basic structure of the functions themselves. 

The key feature enabling a definition of harmonics without separation of variables 
in the Laplacian's eigenfunction equation rests on identifying the Laplacian's symmetry 
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under rotations in d dimensions. Section ^| described these transformations in terms of 
first order infinitesimal rotation operators. We will now construct functions based on 
{i/jj-operators and on the corresponding ladder operators only. Because these functions 
will possess the underlying symmetry of the Laplacian, they are "harmonics" a fortiori. 

The present task amounts to extending the classification of spherical harmonics by 
their label pair (I, m) to multi-dimensional systems, in accordance with procedures to 
construct such harmonics. To this end recall that: 

a) The I label of Y\ m identifies both an eigenvalue of the second-order angular Laplacian 
I 2 and the range < \m\ < I of its second label. 

b) The label m itself is an eigenvalue of the first-order operator l z . 

c) Alternative equivalent sets of harmonics correspond to alternative orientations of 
the z-axis. 

d) Complete sets of spherical harmonics Yi m (6, if) emerge by operating on the 
invariant harmonic Pi(cos8) with the conjugate pair of ladder operators l x ± il y , 
or alternatively operating with the lowering l x — il y alone on the "highest weight" 
harmonic Yu(6,ip). 

Corresponding remarks on hyperspherical harmonics outline here this section's 
development: 

a) We have seen how the multi-dimensional label A performs Z's role in identifying 
eigenvalues of the relevant multi-dimensional angular Laplacian. Alternative 
partitions of A's value will similarly delimit the ranges of m's analogues. 

b) Analogues of m are the integer (or half-integer) eigenvalues mj of the t first- 
order commuting operators {Hi} introduced in section ^|. These sets are viewed 
as components of a vector m in the space subtended by the operators {Hi}. 
They were noted, however, in section p.l| not to be fully compatible with the 
corresponding parameter set {k} that delimits each {|mj|}'s range, equation (|3.3p 's 
label n replacing the eigenvalue m 3 of a different equation. Harmonic eigenvectors 
of the {ifj}'s will be nevertheless identified by convenient sets of m components, 
i.e., by lattice points in the ^-dimensional {Hi} space. 

c) Alternative equivalent sets of hyperspherical harmonics correspond to alternative 
orientations of a vector A in the relevant {Hi} space and to alternative analyses of 
a system's dynamics. Transformations of coordinates and/or of the {Hi} set yield 
equivalent sets of hyperspherical harmonics. 

d) Hermitian-conjugate ladder operators, analogues of physical-space's l x ± il y , will 
emerge as superpositions E a of Hermitian pairs — as anticipated in section — with 
vector labels ex., each label with unit-magnitude components in the {Hi} space, 
being thus represented by diagonal vectors in that space. Linearly independent 
^-dimensional sub-sets of these operators, denoted by {E Va } = {E_ v }, s = 1, . . . , I 
suffice for the present task [7], ||, as well as for encompassing all the far more 
numerous ladder operators by their own appropriate superpositions. 
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A sample set of harmonics emerging in this framework will be displayed at the end of 
section fO| . 

Quantum numbers m^, eigenvalues of Hi operator sub-sets, serve to classify 
"(quasi) -invariants" of multi-particle systems. Such is the J z component of the invariant 
angular momentum J of isolated systems. Two m 8 components pertain to a molecule 
rotating about a symmetry axis of lower inertia in its "body-frame". Three m,'s 
characterize atomic Rydberg states, one of them pertaining to their inner core, one 
to the Rydberg electron and one to their vector sum. A plethora of such m^'s might 
pertain to a turbulent fluid. 

In the absence of quasi-invariants, other than the total «7's, the single label 
m of (|3.7|) may suffice, but a number (< £) of additional m^'s, judiciously chosen 
with reference to the system's structure, helps classifying harmonics Yx^R), and the 
corresponding multi-particle wave functions. 

Geometrical elements of classification also emerge from the nodal patterns of 
harmonics, as noted in section || The Laplacian's separability into coordinates Xi, 
stressed at the outset of section |Q| , affords real solutions of the several resulting one- 
dimensional equations in Xj (with appropriate boundary or periodic conditions) to be 
characterized by n { nodes, yielding altogether £\ rii nodes, a total basically equal to 
the eigenvalue A. Real hyperspherical harmonics with equal A differ then by their A's 
partitions into the relevant n;. 

For the specific purpose of constructing sets of hyperspherical harmonics, sets of 
E± Vs operators, analogues of the physical-space "ladder operators" l x ±il y , complement 
sets of Hi conveniently, much as the l x ± il y complement l z . To this end, the total 
angular momentum's J z = Hi may be complemented by appropriate Hi and thence by 
corresponding ladder operators {E Vs } = {E_ v }. In principle, quantum numbers m s (of 
either sign) of the desired harmonics represent then the number of E± Vs having acted on 
the invariant harmonic C%(R ■ R'), (|3.4j ). In practice, the derivation of hyperspherical 
harmonics for a given A proceeds more appropriately by acting on the analogue of the 
spherical Yu(8, (p) with sequences of lowering operators E^ v , since the "highest weight" 
hyperspherical harmonic is uniquely defined. 

The m s values raised or lowered by E± v operators in this procedure are delimited 
by the relevant "highest weight" eigenvalue A of J z = Hi, viewed as a vector A directed 
along a particular axis and expanded as 



The partition coefficients A s are integers (or half-integers) insofar as both A and the r) s 
are vectors of the {Hi} space with coefficients of integer (or half-integer) magnitude. 
The X s values serve thus as "highest weights", analogues of (|3.3|)'s Zj, for the m s quantum 
numbers. The m s themselves are viewed as components of a "magnetic vector" 




(4.1) 



s 




(4.2) 



s 
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the vectors r) s and hi (the latter pointing in the direction perpendicular to the plane 
of rotation identified by Hi in section 0) being themselves interrelated by integer (or 
half-integer) coefficients. 

An analogue of the spherical harmonics equation defining the range of m, 

(l x ± il y )Y l>±l (e, ip) = (l x ± iy m P,(cos 6) = 0, (4.3) 

is formulated for a hyperspherical harmonic Y\ sTt (R), whose degree is highest (i.e., can 
be raised no further by the ladder operator E v ), in the form 

E v Y XsVa (R) = 0. (4.4) 

When E Vs is cast as a first-order differential operator, (fOj) determines that particular 
"highest weight" harmonic. 

Operating on each of these "highest weight" harmonics with successions of 
"lowering" E_ v operators generates complete sets of hyperspherical harmonics, as 
detailed in the following sections. 

4-1- Operations on the irii parameters 

The mi quantum numbers, eigenvalues of the operators Hi, belong in the hi "space", 
being raised or lowered in value by E a or E± Vb operators external to that space, just as 
the m eigenvalues of l z are shifted by l x ± il y operators with axes orthogonal to z. As 
the ladder operators l x ± il y of physical space are viewed as "eigenvectors" of l z through 
the commutator equations [l z ,l x ± ily] = ±(l x ± ily), sets of ladder operators E a are 
defined as solutions of the commutator-equations' set, 

[H t ,E a ]=a t E a , « 4 = ±lor0; for i = 1, 2, . . .£, (4.5) 

implying that HiE a Ui = (rrij + ai)E a Ui for the eigenvector Ui of Hi with eigenvalue mj. 

Basic elements for solving ( f4.5| ) emerge from the introduction and properties of 
infinitesimal rotation operators ( |2.2|) and ( |2.6| ) in section 0: 

a) Commutators [Hi, Jjjf] vanish unless one of the variables (xj,yk) is common to both 
operators. If one is, the commutator equals a related Jj,^, times a unit-magnitude 
coefficient, much as the three-dimensional commutator of l z with l x or l y does. 

b) Accordingly each E a reduces to a linear combination of a few Jjj! with \ati\ = 1 or 
0; ai = ±1 thus corresponds to raising or lowering mj. The commutator of a pair 
of Hermitian-conjugate ladder operators [E a , E_ a ] is itself Hermitian, specifically 
a multiple of the unit operator. 

c) The directions of a (or r) s ) vectors in the hi space run in that space at equal 
distances between pairs of hi's. 

d) For d even, the resulting operators E a (or E Vs ) shift the values of two mi parameters 
by unity simultaneously, thus preserving the parity of the £\ m«. The same holds for 
most of the ladder operators for odd- dimensional systems, too, with the following 
exception: Odd-dimensional systems include a single operator E Vs (one among a 



Geometry and symmetries of multi-particle systems 



25 



sub-set of i operators E a ) acting on a single m.j, whose contribution violates the 
parity conservation. 



Appendix A outlines a procedure to identify raising and lowering operators. 

These elements, developed originally in [[| |6|], afford a basis for constructing 
complete sets of hyperspherical harmonics by transforming a single prototype harmonic 
with complete sets of E a (or E Vs ) operators. Preferred prototypes have a single non- 
zero value of rrii parameters, typically {mi = A, rrij^i = 0}. An initial E a operator 
will lower mi's value by one unit, raising that of one l^yil from to 1. Alternative 
successions of analogous operators generate thus complete sets of A-degree harmonics, 
all of whose vectors {m} prove compatible with the initial m = {A, 0,0, ..0}, i.e., with 
Y2i \ m i\ — ^ an d (i n the even- dimensional case) (— — (— 1) A - An initial prototype 
harmonic convenient for this purpose reads, in the notation of (|3.3| ), 

cos A a sin A 9^ = ( Xl + iy A * , (4.6) 



R 

with "highest weight" exponents A corresponding to l\ — vn\ — A, l 2 = m 2 = 0, n = 
in Q. 

The allowed range of each among the several m components m 8 is limited by the 
condition \rrii\ < Aj, each Aj being an element of the partition 

A = ^ Ai, A, > 0, -Ai <mi< A;, (4.7) 

i 

which reflects, in turn, specific sequences of the relevant E a applications. Two (or 
more) of the m "weight-vectors" thus generated may coincide, being labelled "multiple 
weights", reached by alternative operator sequences equivalent in this respect. The 
total number of alternative m vectors generated by the present procedure, denoted by 
w in section ^|, includes combined contributions of multiple weights. These w values, 
representing the dimension of the accessible {m} set and of its corresponding set of 
hyperspherical harmonics, are discussed for d- dimensional rotations in f82| , section 10.2] 
and for various other groups in |7|, ||. 

The harmonics' symmetry under rotation reversals, discussed in section [2.1] , implies 
for even-dimensional spaces (where the parity conservation mentioned in item d) above 
holds strictly) the relationship of hyperspherical harmonics 

y\-m{R) = (-l) x Y* m (R). (4.8) 
4-2. Operations on the m s parameters 



Equation (|4.2|) shows parallel expansions of a "weight-vector" m into the eigenvalues 
{mi} and {m s } of the commuting operators {H^} and of the ladder-operator labels 
{t7 s }. These alternative base sets span the same ^-dimensional space with different 
orientations and different metric scales set, respectively, by the Hi eigenvalues nii with 
unit spacing and by the r] s vectors with (in general) two unit-magnitude components 
and resulting squared-magnitudes rj s ■ rj s = 2. Section (O has stressed how ladder 
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operators shift pairs of m; parameters simultaneously, contrasting with E± v 's shift of 
a single m s parameter. The simpler action on the m s thus simplifies the construction 
of hyperspherical harmonics' sets. 

This simplification is partly compensated by a restriction imposed on {m s } sets by 
their equivalence to {mi}, implied by ( |4.2j ) and represented by requiring 

2 HLJk = rs _ qsy (4.9) 

whose non-negative integers r s and q s specify respectively how many times in direct 
succession the operators E_ Va and E v can be applied to m. [For details on this 
equation, see section 10.1 of 0, or sections 13.5 and 15.2 of ||.] 

Two critical elements underlie ( |4.9| ): (i) Equation ( |4.2j ) establishes a linear 
relationship between the parameter sets {mi} and {m s }; (ii) The last relation in ( |4.7p 
restricts the range of each m; sharply. Equation ( |4.9| ) restricts then the ladder operators' 
E± Vs action on ( }4.6| ) or on any hyperspherical harmonic, shifting the relevant m s value 



by one unit, through a sharp selection rule: The resulting value of m s must satisfy (O). 

Within this framework one constructs complete sets of hyperspherical harmonics in 
d dimensions by: 

a) Selecting a coordinate system according to section Q 

b) Constructing an appropriate set of t commuting operators {Hi}; 

c) Constructing corresponding sets of £ ladder operators {E Vs } and {E_ Vg = E^ }, in 
first-order differential form; 

d) Determining a "highest weight" hyperspherical harmonic, characterized by a vector 
A with components Aj in the {Hi} basis and X s in the equivalent {f7 s } basis, by 



solving the relevant ( f4.4| ). [The prototype "highest weight" harmonic ( [4.6| ) pertains 



to A = {A, 0, 0, . . .} in the {Hi} basis.]; 

e) Applying to this "highest weight" harmonic (sequences of) the lowering operators 
E_ Vs , with 1 < s < £, for a total of 2A S times each, to yield the succession of 
harmonics Y\ sms (R), terminating at m s = —X s , remaining however within the 
constraints of ( |4.9|) which often prevents lowering one m s value ahead of other 
m s /'s. 

4-3. Sample derivation of harmonics' sets 



We apply here the procedure just outlined to the three-particle system of section |3T| , 
forming an 5*0(6) geometry, utilizing the same notation and thus taking care of the 
prescription items a) to c), except for selecting a set of rj s vectors. This set, with 
components indicated by a« in (|4.5|), is conventionally M chosen as 



rh~ {1,-1,0}, 77 2 ~ {0,1,-1}, Tfe~ {0,1,1}, (4.10) 
thus implicitly relating each m s number to the mi in (^T). Note how the first of these 



vectors is not orthogonal to the following orthogonal pair, and the set of three is linearly 
independent. 
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We set the value of A in ( f4.6| ) at 2 for simplicity, thus fixing the \ s values in ( |4.1| ) 
at {2, 1, 1}. The relevant Cartesian ladder operators corresponding to the set of {ifj}'s 
r2f) take the form 



E 



(xi ± h/i 
1 



d 



[x 2 =F 12/2J 



5 



<9(x 2 ± iy 2 ) v =f ij/i] 

[jg + j\i =f i (jri - j]?)] 
9 



2 

[zi =F i^2 



5(^2 =F m) 



(x 2 ± ij/ 2 ) 



2 [*^12 + ^22 i i (^12 ^22 )] • 



£7. 



±T ?3 



— 1 



!>1 ± 12: 







-2 



- (ac 2 ± ij/ 2 ) 



d 

d{z 1 ± iz 2 ) 



5 



<9(a?2 =F 12/2) ' d(zi =F Lz 2 ) 

' ^22 i i (^12 + ^22 )] 
Introducing, for the sake of compactness, the elementary operator notations, 



zx 
12 



A + = COS Vi^nr ± 1- 



' <9#j sin 9i difi 

d . „ d 
+ lcot °iTT~ 



1,2: 



we recast ( [4.1 lob e) in polar coordinates, corresponding to {i?j}'s of ( |2.5| ): 



±*7i 



E 



±V2 



E 







sin 9 1 sin #277 tan a sin 9 2 K+ + cot a sin 9\K. 

da ~ 



±itp 2 


cos 


- iii 2) " 


±itp 2 


COS 


+ iL® 



d_ 

'da 



d_ 

'da 



_d_ 

'dfh 



_d_ 

'dth 



'(2) 



'(2) 



(4.11a) 



(4.116) 



(4.11c) 



(4.12 a) 
(4.126) 

(4.13a) 
(4.136) 



(4.13c) 



The prototype harmonic (|4.6| ) with A = 2 reads as cos 2 a sin 9 l e l2ipl ) with A, and \ s 
components {2,0,0} and {2, 1, 1}, respectively. Operators E_ Vs can be applied to this 
expression directly, with full attention to the condition (|4.9| ). This condition excludes at 
the outset acting on this harmonic with either lowering operator E_ V2 or E_ v . i , which 
change mi =2 and mj = 3 from their initial value without affecting m.j = i that retains its 
highest value 2, thus violating the limitations on the mf, only E- n operates on the 
harmonic ( |4.6|) correctly. Sets of harmonics are grouped into "layers" according to the 
number of lowering operators that have been applied to the "highest weight" harmonic 
in the process. Successive steps of lowering the m s quantum numbers are similarly 
restricted, but afford alternative actions of lowering operators, as displayed in table ffl. 
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We actually show the first five layers of harmonics, with the J2 S m s > 0, the harmonics 
with negative values being obtained from these by complex conjugation according to 



Table [I] also illustrates the effect of non-vanishing commutators between ladder 
operators: different ordering of the same ladder operators may yield different harmonics 
with the same label m, i.e., a degenerate eigenvalue, such as m = {0,0,0} in table 
Furthermore, while the procedure demonstrated in table will give the correct number 
of linearly independent harmonics, these may not necessarily be orthogonal (with the 
usual definition of a Hermitian scalar product as an integral over the relevant space: 
(f,d) — J f* #dcc). For a more detailed discussion of the issues of multiple eigenvalues 
and orthogonalization of harmonics in the general case, see e.g. || |84| . 

Finally, recasting {Hi} from (|2.2|) as 

d . , d 



H = -i 



H 



d(xi + iyi 
d 



- (Xi 



d(xi - iy { 
d 



1,2: 



- {Zt - 1Z 2 



(4.14) 



d(zi + iz 2 ) v± * y d(zi-iz 2 ) 
and using ( [4.11 oj -c), we note that manipulating harmonics by means of ladder operators 
is most easily performed in generic Cartesian coordinates, where the harmonics take the 
form of products with generic factors [{xi ±iy i )/R]\ m *\ Q. 



4-4- Harmonics for application to atoms and molecules 

Different sets of hyperspherical harmonics serve to represent different aggregates, at 
different stages of their development. More specifically, each set reflects the structure 
of the relevant "Jacobi tree" introduced in section |3.2| , whose "growth" mirrors the 
integration of the relevant ( |1.2| ), as a function of its hyper-radius R. We deal here 
with particular aspects of harmonics selection, beginning with the remark that, at 
low R values where centrifugal potentials prevail, governing the single particles, the 



set described in section 3.1 may prove adequate, complemented possibly by the set of 



section |3.2| . The following sections deal with two particular aspects of our subject. 



4-4-1 ■ Symmetrized harmonics Particle aggregates include generally sub-sets of 
identical particles: electrons, of course, but also atomic nuclei such as the numerous 
protons of hydrocarbons. As a preliminary to the eventual requirement of anti- 
symmetrizing wave functions under permutation of fermion positions, or analogous 
operations on bosons, it often proves convenient to select at the outset harmonics' sets 
that are invariant under permutation of all identical-particle variables ETl |4"2"f. This 



preliminary operation serves particularly to reduce the dimension of each harmonics' 
set. 

Typically each electron pair may rotate about its centre of mass only with orbital 
momenta equal to an even (odd) multiple of h when in a singlet (triplet) spin state. The 
same holds familiarly for the nuclear states of para- (ortho-)molecular hydrogen. This 



Table 1. Hypersphcrical harmonics for a three-particle system (d = 6) and generalized angular momentum A = 2. The functions are 
constructed by acting on the "highest weight" harmonic (cos a sin #i exp[itpi]) A with different sequences of ladder operators, as indicated in 
the third column. The first two columns list the m quantum numbers in the r/ s and Hi bases, respectively, while the last column identifies 
the harmonics as superpositions of functions (3.3) by listing their corresponding labels {n, l\, mi, I2, m-i). Indentation in the the last two 
columns indicates a continuation line of a single entry 
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restriction reduces the relevant sets' dimensions by exponential factors when applied to 
large sets of identical particles, complemented by enforcing the corresponding symmetry 
between different particle-pairs, a more laborious procedure known as the selection of 
"fractional parentage" [^5| , Chapter 8] , but applied more conveniently at the outset 
of any calculation. 



4-4-2. Expansion of whole state representations Expanding the whole solution ^(R) 
of (|1.2j ) into hyperspherical harmonics may also serve to illustrate the resulting 
representation of a multi-particle state. The set of hyperspherical harmonics must, 
however, be complemented for this purpose by a harmonic function of the hyper-radius 
R, a generalization of the familiar Bessel functions in two dimensions and of their related 
"spherical Bessel" functions in three dimensions. The required harmonic functions of R 
belong once again to the Bessel-function family. 

The spherical Bessel function equation for three dimensions, (10.1.1) of ||83|| , differs 
from the standard Bessel function equation for two dimensions (9.1.1), by: (i) a 
coefficient 2 = d — 1 inserted before its first-derivative term, and (ii) its eigenvalue 
n(n + 1) = (n + 1/2) 2 — 1/4 replacing the standard eigenvalue v 2 . Similarly, the 
hyperspherical Bessel function for d dimensions differs by (i) a coefficient d — l = 3N — 4 
inserted before its first-derivative term, and (ii) its eigenvalue A(A + d — 2) + (d — l)(d — 
3)/4 = (A + (d — 2)/2) 2 — 1/4 replacing v 2 . The resulting Bessel function will thus be 
of integer or fractional order A + id — 2)/2, with a corresponding pre-factor arising from 
the i?( 3Ar_4 )/ 2 coefficient of our equation (jl.2j), as anticipated in |67 . 



5. Hyperspherical expansion of the wave equation 

The preliminary treatment of hyperspherical harmonics in section ^ suffices to formulate 
an expanded version of ( |1 . 2|) . On the left-hand side of this equation, the factors _R( 37V_4 )/ 2 
have been separated out to allow setting a finite initial value of \I/(-R, R) at R — 0. The 
factors separated out take into account the centrifugal — actually, wave-mechanical — 
potential generated in polar coordinates with small values of R by compressing particles 
within short "parallel circles" . We may then standardize wave functions that start at R 
= with unit value in a single hyperspherical channel (Ao, A* ), expanding as R increases 
into alternative channels (A, /x), as indicated by 

A,/x (5.1) 

F\ fj, ,Xfi — > -R A °<5a mo:V' as i? 0. 



Entering the expansion ( |5.1|) in ( |1.2|) , and projecting the result onto the several 
harmonics Yx^R), reduces ( |1.2|) to the system of coupled radial Schrodinger equations 

^Aomo.AmW + E F WvW(fl)l^> = 0, (5.2) 
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with the wave-number matrix 

(AVimiAl.) - (2^- A(A + "- 2)+ J' i - 2)2 - 11/4 ) 
2M(\'n'\Z{R)\\n) 



(5.3) 



i2 

Here Z(R) /R represents the Coulomb potential energy of the interacting particles, 
evaluated at each hyper-radius R. The dots at the end of ( |5.3j ) stand for any additional 
terms of the k 2 matrix corresponding to Hamiltonian terms that represent spin-orbit 
or relativistic corrections not included explicitly. The ability of ( |5.2| ) to include such 
effects — thus far not exploited — may by-pass the current need to treat such terms 
perturbatively rather than directly in the basic equation. 

The system of coupled equations Q5.2|) is formally infinite, owing to the infinite 
range of its parameter A, thus seemingly impractical as noted in section [l|. However, 
circumstances also described in section [l] reduce its size generally to a modest level. 

5.1. Displaying the evolution toward fragmentation 

Our approach to displaying this evolution derives from features of the fragmentation of 
nuclei that are held together by short-range interactions |86||. Beyond this range, r , 



energy eigenf unctions resolve into fragmentation eigenchannels, labelled here by p, each 
of them propagating at r > r in force-free space with unchanged structure, i.e., with 
spherical wave-fronts, r-independent angular distributions f p {6, if) and uniform phases 
4> p (r). The parameter sets {f P {0, ip)} and {<f> p (r)} embody here the effect of all particle 
interactions at r < r . 



The opportunity occurred in [gg] to utilize an analogous parametrization regardless 
of interaction ranges, by embodying the effect of all interactions at the ranges < r < R 
into parameter sets {f p (R; R)} and {<fi p (R)} to be evaluated for successive values of R. 
These sets' dependence on R displays each eigenchannel's evolution as R increases from 
its value (at the system's centre of mass) towards oo, i.e., disregarding all interactions 
at r > R at each step of integration. 

The hyperspherical channel functions ( |5.1| ), ^\ ofJ , (R), of each multi-particle system 
serve here as a basis to calculate the {f p (R;R}} and {(f> p (R)} parameters by casting 
them into superpositions 

$ p (fl) = ^<p(i?)|A /x )* AoMo (J*) = J>(*0|AoA*o> J2 F ^M R ) Y ^( R ) (5-4) 
with initial values of the coefficients 

(p(R)\\olJ, ) -> 1, i.e., (p\ -> (AoMol asR^O. (5.5) 

Requiring each phase (fi p (R) of $ p (i?), and its gradient (dcp/dR), to be uniform 
over each hyper-surface (R = const.) identifies each superposition ( |5.4| ), $ p (i2), as an 
eigenchannel of the propagating Schrodinger equation (|1.2j) at each value of R. It thus 
implies that the angular eigenfunction f p (R;R), representing the aggregate's "shape" 
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at R, evolves "in step" from each hyper-surface to the next. This requirement translates 
into a system of linear homogeneous equations for the coefficients (p(R)\X fj, ) . 



Following the "phase-amplitude" approach of [86], which replaces second order 
wave equations with pairs of first order equations, reference |77|, [78] introduced the 
eigenphase's tan0 p (#) — without previous reference to fl5.4| ) — as an eigenvalue of the "R- 
matrix" 

(\»\K(R)\W) = £ (te) ' *WV(*) 

= ^(A/x|p( J R))tan0 p(/?) (p( J R)|AV). (5.6) 



Reference |79j proceeded then by taking the derivative of ( |5.6|) with respect to R, 



replacing its element (d 2 F/dR 2 ) with its expression (|5.2|) in terms of the k 2 matrix 
(|5.3|), and finally transforming the result to the (p(R)\ basis by means of the coefficients 
(p(R)\Xfj,) and their reciprocals. 

Taking the 1Z matrix as a stepping-stone served thus to resolve the initial second- 
order equation ( |5.2j ) into the set of first-order equations corresponding to the diagonal 
and off-diagonal elements of (|5.6|), respectively, 



dtan0 p(iJ) n ,, , n .,,2, r ,, , n ^ * - 2 



dR 



(p(R)\k 2 (R)\p(R)) t&n 2 cP m , (5.7a) 



d(p{R)\\n) _ ,y sm(j) p{R) {p(R)\k 2 (R)\p'(R))sm(j) pliR) , 
dR sin(</> p - <p pl ) 

with R expressed in units of (dR/d(p p ). This set of equations has been integrated 
numerically in J75| for the prototype example of doubly excited He, with conclusions 
described below. 

5.2. Illustration and discussion 

The numerical integration of ( |5.7a| -fc) for He has been carried out with a device 
that accelerates the convergence of expansion into harmonics (Xfi\, by replacing these 
harmonics with eigenvectors of the k 2 matrix ( p.3| ) at each R. These eigenvectors 
represent "adiabatic" solutions of our problems, carried out earlier in the frame of |T], Q|, 
i.e., disregarding the coupling between the radial and angular variables, R and R. 

The sample results shown in figure 1 should be viewed as interconnecting each 
system's sets of compact (c) and fragmented (/) channels, outlined in |3| and [7^ . 
Quantum mechanical scattering theory represents this connection by "Jost" matrices 
Jf c [[3711 . Each standing- wave eigenfunction of ( |1.2|) , identified by an initial boundary 



condition (c) at its compact limit, is represented asymptotically near its (/) limit by 
sin(kfR)Jf c , or more conveniently in terms of its outgoing and incoming components 

J2 exp(% J R) J+ J2 eM-rt f R)Jj c . (5.8) 
/ / 



Geometry and symmetries of multi-particle systems 



33 




Figure 1. (a) Eigenphases 4> p modulo it, vs the square root of the hyper-radius, 
calculated at the total energy E — 57 eV above the ground state, just below the 
"(2s 2 ) 1 <S"' resonant state in helium, (b) shows the same set of eigenphases, calculated 
at a total energy E — 58 eV above the ground state. 



The matrices Jf c are complex-conjugate for "open" channels /, i.e., when energy 
suffices to achieve the / limit; for energetically "closed" channels the wave-number kf 
is imaginary, whereby one component converges to zero at discrete eigenvalues of the 
energy E and the other diverges as R — > oo. The Jost matrices serve then to construct 
scattering matrices 

s ff = E J h (j=) f > ( 5 - 9 ) 

as detailed in || In this frame one views each amplitude (p(R)\X fj, Q ) of (|5.4j) , 
evaluated at a finite range R, as a partial construction of the Jost matrix element J/ c 
with / representing the limit of (p(-R)l as R ~ > 00 • Plots of the several bra symbols 
(p(R)\ achieve our objective of displaying the system's evolution from R = toward oo. 

Figures |l](a) and (b) plot eigenphases 4> P {R) vs. y/R, modulo rr, at energies 
straddling the (2s) 2 1 S resonance of He near 58 eV, for a number of (p(R)\Xo/j, ) pairs. 
Each line's slope mirrors the rate of increase of <f> p (R), i.e. (loosely), the rate of expansion 
in R of the corresponding eigenfunction. [The first eigenphase's slope reflects the rapid 
motion of an electron ionized with approximately 33 eV kinetic energy; successive curves 
reflect the increasingly slower development of two-electron excitations in successively 
higher modes.] In the lower and upper ranges of the ordinate (f) p , corresponding to low 
values of sm[4> p (R)], pairs of curves appear to cross with minimal disturbance, owing 
to unresolved values of the coupling coefficient on the right of ( |5.7t|) , in spite of the 
singularity arising from the vanishing of its denominator at each crossing. Major effects 
of crossings emerge instead at middle ranges of <p p (R), where pairs of curves appear to 
repel each other experiencing major deflexions. 

The localization of such major "avoided crossings" within limited ranges of 
ordinates and abscissas, confirmed by analogous more dense plots, verifies our 
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expectation that eigenfunctions propagate smoothly outside the limited parameter 
ranges where the coupling coefficients of ( |5.76Q diverge. The apparent "repulsion" 
of curves at avoided crossings is a familiar aspect of the Landau-Zener phenomenon 
f88| , |89| , p0|1 occurring at points of "degeneracy" where two alternative values of a 
parameter coincide, as the phases (0 P , </y) do here (modulo 7r). 

An additional major manifestation of avoided crossings, namely, the hybridization 
of pairs ((p|A At ), (//|Aq/Xq)) has, however, not been included in the calculations 
underlying figure [l]. This hybridization might result by fitting parameters of the observed 
phenomena — slopes, slope differences of the curves and their closest approaches — to 
the corresponding elements of the Landau-Zener theory. The relevance of the present 
(p|A £t ) coefficients to the Jost matrices of interest remains fragmentary pending 
further developments of Nakamura's recent analytical improvements on the Landau- 
Zener theory g|. 



5. 3. Qualitative effects of Coulomb interactions 



Equation (|5.3|) 's matrix (X'fj,'\Z(R) |A/x) consists of N(N — l)/2 terms for an aggregate 
of N charged particles. These terms, included sequentially in computer programs, serve 
to solve (|5.7q|) and (|5.76| ) numerically, yet warrant analysis aimed at visualizing their 
qualitative action. To this end one may resolve the Z-matrix' action on the wave function 
\1> into its several aspects. 

Sub-sets of terms acting between particles with equal (opposite) charges push 
corresponding particle sub-sets apart (together). Within this scope one notes that: 

a) Expressing the distance between each atomic nucleus and another particle in terms 
of their common mass-weighted hyper-radius R scales up its charge by its mass' 
square root, thereby boosting its interactions' strength and thus favouring molecular 
dissociation as compared to ionization by electron ejection. 

b) The force acting on each particle-pair depends on the pair's orientation. 
Combining the Coulomb forces between various particle-pairs involves thus 
extensive geometrical transformations of the relevant position coordinates. 

Coulomb interactions between sub-sets of identical particles may be conveniently 
grouped, particularly so following the symmetrization of relevant position coordinates 
outlined in section |4.4| and in [Appendix D|. It would then become possible to combine 



Coulomb terms pertaining to such sub-sets to yield expressions of their electrical 
multipole moments and of their corresponding multipole fields. A semi-macroscopic 
view of each aggregate's mechanics should thus emerge, in terms of collective variables. 

Within the context of initial combination of molecular nuclei with closed-shell 
electrons, one could then envisage treating all molecular valence electrons as forming 
a gas whose electric multipole moments are inflated by electronic mutual repulsions, 
yet contained by the attractive multipole fields of nuclei and closed-shell cores. The 
proton-nuclei of hydrocarbon molecules would be similarly viewed. This attractive 
containment would perform a two-fold action: Holding the molecule together as a unit 
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and simultaneously smoothing out the distribution of opposite charges throughout its 
volume. 

Appendix A. Construction of ladder operators 

The generic infinitesimal operators ( |2.6| ) raise and lower the eigenvalues, rrii and rnj, of 
operator pairs (Hi, Hj) by mapping a single harmonic onto a superposition of harmonics 
with different m^'s and m/s. We construct here linear combinations of those operators 
that raise or lower by unity these eigenvalues in a definite way, according to 



Casting for this purpose (H i} Hj) in the form of the first two expressions in (|2.2|) , 
i.e., as (Hi = J??, Hj = Jj?), restricts the solutions E a of fl4.5|) to combinations of the 
four components 

jxx jyy jxy jyx ( A 1\ 

°ij J ° ij J ° ij J u ij ■ K^- 1 ) 

Identification of the proper linear combinations proceeds through analysis of their 
commutator relations with the H^s, namely, 

[Hi, m = • i-C7- Wi, m = ■ l, u- ( A - 2 «) 

[Hi, J!-} = -iJg, [Hj, JJ] = -iJg, (A.26) 

[Hi, J??] = -iJ%, II r J??] = (A.2c) 

[Hi, J%\ = [H S , Jfj] = I U'f. (A.2d) 

Combinations of the operators (|A.1[) symmetric and antisymmetric in their (x, y) 



variables, analogues of (|2.7fj )'s symmetry in (i,j) indices, yield 

E ±c ^ +) = (J*f - Jf?) ± + Jf), (A.3o) 

eJL-) = (J--+J1!)T - J?*), (A.36) 

satisfying the desired equations 

[H, E ± J + -,} = ±E ±a(+ - h [H 3 , E (+-)} = TE ±a(+ - h (AAb) 



with 



og +) = (+l,+l), (A.5a) 
ag- ) = (+l,-l), (A.56) 



c4 ++) = a^ +_) = 0, for k ^ i,j. (A.5c) 

Note that the two operator sets {\(Jf 3 x T Jfj), \(J% ± J f/)> U H j ± 
commute exactly like {l x ,l y ,l z } i n three dimensions. The occurrence of (symmetric 
or antisymmetric) pairs of rotation operators in the role of both l x and l y , respectively, 
reflects the feature of the pairwise change of m, and rrij. The symmetry of these operator 
pairs under interchange of the first and second coordinates Xi <-> yi and Xj <-» yj of the 
relevant Hi and Hj extends the present construction to ladder operators that change 
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the m quantum numbers of H operators of the third type in ( |2.2| ), involving two z- 
coordinates. 

For odd dimensions, t ladder operator pairs pertaining to coordinate pairs (i, 0) 
change only one of the m^s. If again X{,yi make up the operator Hi pertaining to 
rrii, the two ladder operators acting on mi are complex linear combinations of the two 
infinitesimal rotation operators which involve either of Xj or y^, together with the single 
unpaired coordinate [denoted x in section e.g. the z coordinate in the familiar case 
of SO (3)]. In this case, the ladder operators are completely analogous to the ones in 
SO (3). 



Appendix B. Transformation between Jacobi trees 

As an example of transformation between Jacobi trees, consider the two trees 




The mass-weighted relative coordinates for tree A result from independent-particle 
coordinates by the transformation 



£i = a/tt 7-H r 2 — fi) B.la) 

Sl V Mi + M 2 V 1 11 v ' 

A I .U,.U, 



M 3 + M 4 



(Mi + M 2 ) (M 3 + M 4 ) / M 3 r 3 + M 4 r 4 MjTi + M 2 r 2 \ 

M 1 + M 2 + M 3 + M 4 V M 3 + M 4 Mi + M 2 / I • C J 



A = / (Mi + M 2 + M 3 + M 4 )M 5 / _ Min + M 2 r 2 + M 3 r 3 + M 4 r 4 \ 
44 y Mi + M 2 + M 3 + M 4 + M 5 V 5 Mi + M 2 + M 3 + M 4 / ' 1 ' j 

i.e. by first connecting particles 1 and 2, then 3 and 4, then the complex {12} to 
the complex {34}, and finally the complex {1234} to 5. The fifth Jacobi coordinate 
represents the centre of mass which remains the same for all trees consisting of the same 
particles, and is hence irrelevant to transformations of the four relative coordinates. 
The transformation from A to D resolves into three elementary "transplantations" : 
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12345 12345 12345 12345 

Tab \ \ y / Tbc \ V / / Tqd 







The transformation Tab from A to B affects only the first and the third among the mass- 
weighted Jacobi coordinates, since it transplants branch 2 from the complex {12} to the 
complex {234}. It is therefore represented by applying to the four-vector (£ x , £ 2 , £3 , £ 4 ) 
the matrix 

sin (pAB \ 


1 



/ COS (pAB 





- AB 



\ 



sin (pAB 





1 





COS (pAB 





(B.2o) 



/ 



i.e. by a "kinematic rotation" through an angle 



'AB 



tan 



' M 2 (M 1 + M 2 + M 3 + M 4 ) 



(B.26) 



Mx(M 3 + M 4 ) 

according to (|2.13|). Note the general structure in the mass coefficients: Transplanting 



branch q from the complex {pq} to the complex {qr} corresponds to a rotation by 
= tan -1 a/ 'M q (M p + M q + M r )/M p M r (in the first quadrant, i.e. with positive signs 
for both cosine and sine). The transformation Tbc transplants branch 3 from {34} to 
{23}, and Tqd transplants the complex branch {234} from {1234} to {2345}, with the 
respective transformation matrices 



T, 



BC 



\ 



COS <pBC 

sin 4>bc 





- sin (pBc 
cos 4>bc 







\ 



1 
1 



T< 



CD 



( 1 

10 

cos (pcD — sin^CT) 

\0 sm(f) CD cos (pcD j 



(B.3o) 



and rotation angles 



?CD 



tan 



tan" 



1 M 3 (M 2 + M 3 + M A ) 



M 2 M A 

\M 2 + M 3 + M 4 )(Mi + M 2 + M 3 + M 4 + M 5 ) 
MiM 5 



(B.36) 
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The complete transformation for this sequence is represented by Tad 
with 



Tcd ' Tpm • T 



BC " J-AB 



T 



AD 



(B.4) 



-sm.<j)BC sin (pAB cos (pBC \ 

cos 4>bc - sin 4>ab sin 4> B c 
cos 4>ab cos 4>cd - sin 4> C d 

COS 0y4iJ sin (pCD COS (f) C D ) 

Inserting the explicit expressions for the angles <\> verifies that the coordinates 1 < 
i < 4, indeed describe the relative coordinates of tree D in terms of independent particles 
coordinates, namely 



/ COS (pAB COS <pBC 

cos 0ab sin BC 
sin <pAB cos 0CD 
y sin 4>ab sin CD 



1) 



(M 2 + M 3 )M 4 f M 2 r 2 + M 3 r 3 _ 



M 2 + M 3 + M 4 V M 2 + M 



r 4 



/ M 2 M 3 
M 2 + M 3 



(r- 2 - 



(M 2 + M 3 + M 4 )M 5 
M 2 + M 3 + M 4 + M 5 



M 2 r 2 + M 3 r 3 + M 4 r 4 
M 2 + M 3 + M 4 



r 5 



M 1 (M 2 + M 3 + M 4 + M5) / M 2 r 2 + M 3 r 3 + M 4 r 4 + M 5 r 5 
Mi + M 2 + M 3 + M 4 + M 5 V M 2 + M 3 + M 4 + M 5 



ri 



(B.5a) 
(B.56) 
(B.5c) 

(B.5d) 



Appendix C. Finite Transformations of Hyperspherical Harmonics. 

Finite transformations of multi-dimensional harmonics (or operators) correspond to the 
infinitesimal ones considered in section |] just as the prototype transformation (|2.3|) 
(pertaining to a physical-space rotation) corresponds to the infinitesimal Q2.1Q . This 
correspondence holds generally, since all transformations relevant to this paper resolve 
into products of two-dimensional rotations, as stressed repeatedly in the text. 

Section (| has identified hyperspherical harmonics in the frame of a representation 
based on a vector A in the ^-dimensional space of maximally commuting operator sets 
{Hi}. Within this scope, we might deal here just with ^-dimensional rotations of A. 
This space itself is, however, subject to rotations of the {Hi} induced by the 3(N — 1)- 
dimensional coordinate rotations considered in section |2| for TV-particle aggregates. 
Generic infinitesimal operators on such spaces were indicated in (|2.6|) by J^J 1 , whose 
label xy refers to a pair of coordinate axes, whereas ij refers to a pair of particles. 

Viewing physical-space rotations, identified by three Euler angles, as our model, 
recall how two of these angles pertain to rotations about (hence parallel to an 

xy-plane) and the third one to a shift of z's orientation to a new direction z', usually 
understood to lie on the previous xz-plane. Whereas rotations by an angle <p about 
z simply multiply eigenvectors of l z with eigenvalue m by e im</3 , rotations by an angle 
9 in the xz-plane transform it into a superposition of eigenvectors whose eigenvalues 
m! result by transforming the initial m with the Wigner matrix d^, m (6). [The index I 
stands here for the largest value ("highest weight") of m.) 
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Correspondingly, in a multi-dimensional space, we consider two distinct classes of 
two-dimensional rotations: (i) Rotations about one of the {Hi} operators' symmetry 
axes (i.e., in a plane orthogonal to that axis) which multiply a harmonic eigenvector of 
Hi with eigenvalue \rrii\ < Aj by e iniilfii ; and (ii) Orientation changes of an H^s own axis, 
within a specified plane through that axis, yielding a superposition of harmonics with 
Wigner coefficients D J , whose subscripts differ only by replacing their component 
with m\. 

The multi-dimensional framework deals with transforming from one coordinate 
basis, including its operator set {Hi}, to a new basis with its operator set {H'A, each 
set of indices {i = 1,2, . . .} and {j = 1,2,...} being ordered. This framework affords 
articulating generic transformations through sequences of two-dimensional rotations 
by Euler angles {f a } and {0 U }: The initial ip a =i equals the angle between the ifj's 
zero-azimuth and the plane of the pair {Hi=i, iifj =1 }'s axes. The initial Q v =\ equals 
similarly the angle between the axes of H i= i and H'- =1 . The next ip a =2 shifts the plane 
of the H i= i, Hj =1 axes to H'j =1 's zero-azimuth. Corresponding angles {(p a =3, 9 U =2, ^Pa=i} 
pertain to the operator pair {H i= 2,H'- =2 }, a procedure to continue recursively. 

Insofar as the {Hi} operators are anchored to their coordinate systems, their two 
classes of rotations drag their coordinate axes along. Altogether, transformations of 
harmonics indices by coordinate rotations are thus seen to resolve into three elements: 
(i) Rotation of the "representation vector" A in the ^-dimensional space of the {Hi} 
set; (ii) Rotation of the {-ffjj's themselves, described above; (iii) Further rotations of 
coordinates with respect to the {Hi}. 



Example: Kinematic rotation of harmonics 

We transform here hyperspherical harmonics of two Jacobi vectors representing a three- 
particle system. Equations ( |2.12| ) and ( |2.13| ) have described the transformation of a 
Jacobi vector pair from tree A to tree B 



A: I ^V S(r '" r2) (C.la) 

' {A (M 1 +M 2 )M 3 , M iri +M 2 r 2 \ ' 

<S2 V M1+M2+M3 V M1+M2 3 ^ 



B: { Cf_ vS^-^ (ai6) 

1 »B _ hh{M 2 +M s ) / _ M 2 r 2 +M z r 3 \ 
S2 y M1+M2+A/3 \' 1 M2+M3 n 

kinematic rotation by an ang le /3 = tan" 1 y J{M 1 + M 2 + M 3 )M 2 /(M 1 M 3 ), 

ttf , t%} = {a»/9tf - sin/3^, sm(3£? + (C.2) 

Since each vector £ has three spatial components, we deal here with a six-dimensional 
coordinate transformation with components 

{tf, } = {(#)«, (tf )y, (tf )„ (tf )«, (tf )y, (tf ),} = {*!, yi, «1, X 2 , 1/2, ^2}, (C.3) 

and similarly for £f . The x-component of £f results from a rotation by the angle /3 in 
the Xix 2 -plane, the x-component of £f by a rotation through — (3 in the (oriented) x 2 Xi- 
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plane, and likewise for their y- and components. The corresponding transformation 
of three-particle harmonics |A, fj, A ), as represented, e.g., by (|3.3|), is indicated according 
to section |3| by 

\\p B ) = T ab (P)\^a) = E \ X »'a)D x ^ a {P). (C.4) 

Dealing here with Cartesian coordinate rotations of the £ vectors, at variance with 
the preceding polar coordinate description, our present operator Tab factors into three 
separate (commuting) transformations of harmonics corresponding to rotations by (3 in 
each of the three coordinate planes X1X2, 2/12/2, and z\z 2 

T AB ((3) = exp(ipj%) exp(i/?J?f ) exp(i/3Jff ). (C.5) 

Convenient expressions of the three infinitesimal operators in this expression of Tab 
appear in earlier parts of this paper: (i) Equation (|2.2|) identifies Jff as the operator 
H 3 for our three-particle system, yielding immediately 

exp(i/3Jg)\\iJt A ) = exp(im 3 /3)|A/i A ); (C.6) 

(ii) Equation ( |A.3£| ) identifies the combination Jff + Jf| as sum of the operator E (+-) 
and its reciprocal E (+-), both acting on the eigenvalues m\ and mi of the \H\,Hi\ 

pair, the symbol ck^ - '* meaning "raising m\ and lowering 7712". The combination 
E (+-) + E <+-) of two reciprocal (Hermitian conjugate) operators is itself Hermitian. 

a 12 Q 12 

The resulting matrix elements (Xfjf A \ exp(i/3{E {+-)-\-E_ (+-)}) lAua) amount to Wigner 
^m'm elements, with parameters specified by the following observations. Projecting 
Ha an d V>'a on ^° a V2 generalizes the lower indices in a^, m to the present higher- 
dimensional setting: 

m= ^"V ^-^), (C.7) 

a l2 ' a 12 



and similarly m! = (m'j — m' 2 )/2. Section |4.2| discussed in item e) a succession, or 
chain, of harmonics labelled here by \x A + raQ: i2 _ ' ) 5 with total length 2A S , restricted 
however by (|4.9| ). The "multipole order" represented by the upper parameter in d^, m 
corresponds here to half the length of this chain of harmonics. Finally, the generalization 
of the "triangular condition" familiar from three dimensions now requires \\fJ> A ) an d 
|A/x^) to lie on the same chain of harmonics |A, fx A + nOL 12" V Analogy to the relation 
lx = h{l+ + 1-) suggests treating exp[i/3(E (+_ > + E _ (+-))] like a rotation about the 
x-axis in three dimensions. With the standard definition of Euler-angle rotations, the 
corresponding matrix element picks up an additional phase factor exp[i(m' — m)7r/2] 
resulting from rotating the y- onto the x-axis and back. Note also that the sum of 
reciprocal operators in our matrix element lacks the factor 1/2, thereby effectively 
multiplying the angle (3 in the <i-symbol's argument by a factor 2. The procedure, 
outlined here for a three-particle example, extends similarly to larger aggregates. 
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Appendix D. Outline of procedures for treating large sets of identical 
particles 

Constructing wave functions of a few (3-4) electrons, with the required anti- 
symmetrization, is rather familiar, being extended to atomic shell- filling, e.g., in 



Chapter 8 of |82| . Its extension to much larger sets remains problematic. Deceptively 
simple considerations, to be presented below, indicate that this extension may actually 
proceed along the same lines, essentially because successive steps prove independent 
of one another. Whether these considerations constitute more than just a "solution 
in principle" for this fundamental problem of quantum many-body theory remains 
to be seen. The inevitable exponential proliferation of operations to be carried 
out hampers their actual implementation for all but the smallest sets of particles. 
Significant simplifications occur particularly for a system's ground state configuration 
owing to (in general) higher symmetry in this state. Even for this case, however, the 
anti-symmetrization of sets with more than three identical particles still presents a 
formidable task. The most promising approach to this problem's systematic solution 
and its implementation is currently being developed by Barnea and Novoselsky J|!| 
employing the concepts discussed in this Topical Review, namely, Jacobi coordinates 
and hyperspherical harmonics. 



Coordinates and their symmetrization 

Identity of particles implies that permutation of any pair of them leaves any function 
of the pair unchanged. Artificial labelling of such particles, by indices i = 1,2,3,..., 
appears nevertheless generally desirable for purposes of "book-keeping". It is then 
necessary to "symmetrize" any function of particle positions, /(r^, r 2 , . . .) to ensure 
its invariance under permutation of each pair of indices {i, j), by superposing sets of 
such functions differing by the whole set of relevant permutations. Spin and position 
coordinates should be combined, of course, in this construction that proves increasingly 
laborious with increasing number of particles. 



Classification by seniority 

The "seniority" label ( 'V ) of an atomic state indicates the number of its particle pairs 
characterized as l S and thus isolated from its remaining particles. This characterization 
means invariance under rotation of space coordinates of both spin and position variables. 
The spin part of this label applies equally regardless of the total number of particles 
in the system. The "5"' label implies a spherically symmetric matching of the pair's 
angular distribution, whose extension to multi-dimensional systems needs elaboration. 
Considering that this symmetry is attained by combining a pair of orbitals even and 
odd, respectively, under reflection through a plane — thus stretching in orthogonal 
directions — we suggest achieving the corresponding invariance in higher dimensions by 
combining pairs of hyperspherical harmonics pertaining to Jacobi trees constructed by 
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selecting orthogonal space directions at each step-wise addition of one particle. 



Introduction of "triple tensors" 



This operation, introduced by Judd in the 1960's and described in |82|, pp 209ff], yields a 
systematic classification of the shell-filling process for the electrons of each atomic shell. 
It rests on elementary applications of coordinate-rotation transformation pairs and of 
their "reduction" , which appear as equally serviceable regardless of their dimensionality. 



Separation of particle sub-sets 

This operation, familiar in atomic systems and leading there to the "fractional 
parentage" procedure, has a major role in multi-particle settings where, typically, 
electron sub-sets perform varied functions forming "closed shells" (or sub-shells) of 
different atoms as well as chemical bonds, preserving the relevant symmetry and 
coherences. The structure and flexibility of the Jacobi trees corresponding to alternative 
hyperspherical harmonics appear well suited to extension to multi-particle settings, with 
appropriate development of recursion techniques. 



References 



[1] Macek J H 1967 Phys. Rev. 160 170 

[2] Macek J H 1968 J. Phys. B: At. Mol. Phys. 1 831 

[3] Fano U 1981 Phys. Rev. A 24 2402 

[4] Fano U 1983 Phys. Rev. A 27 1208 

[5] Cartan E 1894 Sur la structure des groupes de transformations finis et continus (Thesis) [Reprint 

1952 in: Qiuvres completes (Paris: Vuibert) Part I, Vol I, p 127ff] 
[6] Racah G 1965 Group Theory and Spectroscopy, in: Ergebnisse der exakten Naturwissenschaften 

vol 37 (Berlin: Springer) p 28 
[7] Wybourne B G 1964 Classical Groups for Physicists (New York: Wiley) 
[8] Cornwell J F 1984 Group Theory in Physics (San Diego: Academic), Vol II, Chapters 13-16 
[9] Morse P M and Feshbach H 1953 Methods of Theoretical Physics (New York: McGraw Hill) p 1730 

[10] Fock V 1958 Norske Vidensk. Selsk. Forh. 31 138 

[11] Wannier G H 1953 Phys. Rev. 90 817 

[12] See, e.g., Smith F 1960 Phys. Rev. 120 1058 

[13] Tsipis C A, Popov V S, Herschbach D R and Avery J S (eds.) 1996 New Methods in Quantum 

Theory (Dordrecht: Kluwer) 
[14] Delves L M 1959 Nucl. Phys. 9 391 
[15] Delves L M 1962 Nucl. Phys. 20 268 

[16] Cobis A, Fcdorov D V and Jensen A S 1998 Phys. Lett. B 424 1 
[17] Kievsky A 1998 Nucl. Phys. A 631 669c 

[18] Gorbatov A M and Khazov A Yu 1997 Phys. At. Nucl. 60 1979 

[19] Santopinto E, Iachello F and Giannini M M 1998 Eur. Phys. J. A 1 307 

[20] Smirnov Yu F and Shitikova K V 1978 Sou. J. Part. Nucl. 8 344 

[21] Ballot J L and Fabre de la Ripelle M 1980 Ann. Phys., NY 127 62 

[22] Zhukov M V, Danilin B V, Fedorov D V, Bang J M, Thompson I J and Vaagen J S 1993 Phys. 
Rep. 231 151 

[23] Kuppermann A and Kaye J A 1980 Chem. Phys. Lett. 74 257 



Geometry and symmetries of multi-particle systems 



[24] Kaye J A and Kuppermann A 1981 Chem. Phys. Lett. 77 573 

[25] Kaye J A and Kuppermann A 1981 Chem. Phys. Lett. 78 546 

[26] Pack R T and Parker G A 1987 J. Chem. Phys. 87 3888 

[27] Linderberg J, Padkjaer S B, Ohrn Y and Vessal B 1989 J. Chem. Phys. 90 6254 

[28] Launay J M and Le Dourneuf M 1989 Chem. Phys. Lett. 163 178 

[29] Launay J M and Le Dourneuf M 1990 Chem. Phys. Lett. 169 473 

[30] Lepetit B and Launay J M 1991 J. Chem. Phys. 95 5159 

[31] Honvault P and Launay J M 1998 Chem. Phys. Lett. 287 270 

[32] Groenenboom G C 1998 J. Chem. Phys. 108 5677 

[33] Yu H G 1997 Chem. Phys. Lett. 281 312 

[34] Wu Y M, Lepetit B and Kuppermann A 1991 Chem. Phys. Lett. 186 319 

[35] Adhikari S and Billing G D 1998 Chem. Phys. Lett. 284 31 

[36] See, e.g., Aquilanti V, Bonnet L and Cavalli S 1996 Mol. Phys. 89 1 

[37] Aquilanti V and Avery J 1997 Chem. Phys. Lett. 267 1 



[38] Aquilanti V, Cavalli S, Coletti C, De Fazio D and Grossi G 1996 in: fl3j p 233fT] 

[39] Aquilanti V, Cavalli S and Coletti C 1998 Phys. Rev. Lett. 80 3209 

[40] Novoselsky A, Katriel J and Gilmore R 1988 J. Math. Phys. 29 1368 

[41] Novoselsky A and Katriel J 1994 Phys. Rev. A 49 833 

[42] Barnea N and Novoselsky A 1998 Phys. Rev. A 57 48 

[43] Ovcharenko V I 1976 Sov. J. Nucl. Phys. 24 483 

[44] Lin C D 1981 Phys. Rev. A 23 1585 

[45] Lin C D 1982 Phys. Rev. A 25 1535 

[46] Lin C D 1984 Phys. Rev. A 29 1019 

[47] Chen Z and Lin C D 1990 Phys. Rev. A 42 18 

[48] Lin C D 1982 Phys. Rev. A 25 76 

[49] Lin C D 1982 Phys. Rev. A 26 2305 

[50] Lin C D 1983 Phys. Rev. A 27 22 

[51] Lin C D and Macek J H 1984 Phys. Rev. A 29 2317 

[52] Feagin J M and Briggs J S 1988 Phys. Rev. A 37 4599 

[53] Watanabe S and Lin C D 1986 Phys. Rev. A 34 823 

[54] Tolstikhin O I, Watanabe S and Matsuzawa M 1995 Phys. Rev. Lett. 74 3573 

[55] Christensen-Dalsgaard B L 1984 Phys. Rev. A 29 2242 

[56] Tang J, Watanabe S and Matsuzawa M 1992 Phys. Rev. A 46 2437 

[57] Sadeghpour H R and Greene C H Phys. Rev. Lett. 65 313 

[58] Sadeghpour H R Phys. Rev. A 43 5821 

[59] Qiu Y, Tang J-Z, Burdorfer J and Wang J 1998 Phys. Rev. A 57 R1489 

[60] Greene C H 1981 Phys. Rev. A 23 661 

[61] Watanabe S 1982 Phys. Rev. A 25 2074 

[62] Clark C W and Greene C H 1980 Phys. Rev. A 21 1786 

[63] Bao C G, Yang X and Lin C D 1996 Phys. Rev. Lett. 76 3096 

[64] Yang X, Bao C G and Lin C D 1996 Phys. Rev. A 53 3934 

[65] Morishita T, Lin C D and Bao C G 1998 Phys. Rev. Lett. 80 464 

[66] Morishita T and Lin C D 1998 J. Phys. B: At. Mol. Opt. Phys. 31 L209 

[67] Cavagnero M 1984 Phys. Rev. A 30 1169 

[68] Cavagnero M 1986 Phys. Rev. A 33 2877 

[69] Cavagnero M 1987 Phys. Rev. A 36 523 

[70] Braun M, Schweizer W and Elster H 1998 Phys. Rev. A 57 3739 

[71] Braun M, Schweizer W and Elster H 1998 Nucl. Phys. A 631 644c 

[72] Aquilanti V and Avery J 1997 Chem. Phys. Lett. 267 1 

[73] Fano U 1983 Rep. Prog. Phys. 46 97 

[74] Lin C D 1986 Adv. At. Mol. Phys. 22 77 



Geometry and symmetries of multi-particle systems 



44 



[75] Lin C D 1995 Phys. Rep. 257 1 

[76] Chattopadhyay R and Das T K 1997 Phys. Rev. A 56 1281 

[77] Fano U and Sidky E 1992 Phys. Rev. A 47 4776 

[78] Bohn J L 1994 Phys. Rev. A 49 3761 

[79] Bohn J L and Fano U 1996 Phys. Rev. A 53 4014 

[80] Avery J 1989 Hyper spherical Harmonics (Dordrecht: Kluver) 

[81] Erdelyi A, Magnus W, Oberhettinger F and Tricomi F G 1953 Higher transcendental functions 

(New York: McGraw-Hill) Vol II, Chapter 11, p 232ff 
[82] Fano U and Rau A R P 1996 Symmetries in Quantum Physics (San Diego: Academic) 
[83] Abramowitz M and Stegun I A 1965 Handbook of Mathematical Functions (New York: Dover) 

Chapter 22 [Ninth printing 1972] 
[84] Heim T A and Green D 1999 J. Math. Phys. at press 

[85] Cowan R D 1981 Theory of Atomic Structure and Spectra (Berkeley: University of California 
Press) Chapter 9 

[86] Calogero F 1967 Variable Phase Approach to Potential Scattering (New York: Academic) Chapter 
19 

[87] Newton R G 1966 Scattering Theory of Waves and Particles (New York: McGraw Hill) 

[88] Landau L D 1932 Phys. Zts. Sowjet. 2 46 

[89] Zener C 1932 Proc. R. Soc. A 137 696 

[90] Stiickelberg E C G 1932 Helv. Phys. Acta 5 369 

[91] Zhu C and Nakamura H 1992 J. Math. Phys. 33 2697 

[92] Zhu C and Nakamura H 1995 J. Chem. Phys. 102 7448 



